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SECTION  I 
INTRODUCTION 


The  application  of  computers  in  research  programs  is  rapidly 
increasing  due  to  the  inherent  advantages  possible  in  the  use  of 
these  machines.  Experiments  that  were  impractical  or  impossible 
with  human  operators  are  now  feasible,  mathematical  models  can  be 
constructed  and  readily  evaluated,  sophisticated  data  reduction 
involving  millions  of  calculations  can  be  performed  in  seconds  to 
enhance  data  analysis,  computer  graphics  are  providing  an  inter- 
active dialogue  with  the  computer  to  evaluate  model  building.  All 
of  these  capabilities  are  being  exploited  to  a greater  and  greater 
degree  as  one's  awareness  of  the  computer's  capability  increases 
and  the  relative  time  and  cost  per  calculation  continues  to  decrease 
The  Air  Force  Materials  Laboratory  use  of  computer  has  also  been 
increasing.  While  the  majority  of  the  past  computer  experience  of 
the  Materials  Laboratory  has  been  with  batch  jobs  to  a large  com- 
puter such  as  the  CDC  6600  and  high  cost,  specifically  designed 
minicomputer  systems,  the  introduction  of  the  microprocessor  has 
provided  a flexible,  low  cost  solution  to  obtaining  maximxam  effi- 
ciency in  in-house  and  on-site  research  programs. 

\/This  report  summarizes  the  work  performed  in  pi'oviding  scien- 
tific and  engineering  computer  support  to  various  projects  within 
AFML  during  this  time  period.  This  work  includes  exploitation  of 
all  computers  available  to  AFML  researchers.^?  While  the  results  of 
this  work  have  been  or  are  in  the  process  of  being  published  in 
either  technical  reports  or  scientific  journals,  completion  of 
some  projects  and  publication  of  the  results  will  occur  on  a suc- 
ceeding contract. 

^^This  report  is  divided  into  major  topics:  (1)  Automated 
Experiments  and  Data  Acquisition,  (2)  Computer  Programs  for  Re- 
entry and  Laser  Applications  of  Ablation/Thermal  Effects,  (3)  Com- 
puter Analysis  of  Microstructures,  and  (4)  General  User  Programs 
and  Education, 
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SECTION  11 


AUTOMATED  EXPERIMENTS  AND  DATA  ACQUISITION 


2.1  AFML  DIGITAL  COMPUTER  HIERARCHY 

It  is  the  object  of  the  Computer  Activities  Office  of  the  AF 
Materials  Laboratory  to  organize  a digital  computer  hierarchy  so 
that  with  the  automation  of  laboratory  experiments,  data  from  the 
experiments  can  be  transmitted  to  other  digital  computers  within 
ASD  for  data  reduction.  A block  diagram  of  this  hierarchy  is  shown 
in  Figure  1.  This  communications  requires  the  use  of  telephone 
communicaticns  equipment  such  as  modems,  acoustic  couplers,  audio 
telephone  lines,  and  hard-wired  telephone  lines.  For  the  Materials 
Laboratory,  the  central  point  of  this  hierarchy  is  the  SEL  86  located 
in  Building  652.  Data  from  the  experiments  is  gathered  by  the  mini 
and  microcomputers  and  transmitted  to  the  SEL  via  either  low  speed 
(110,  300  Baud)  or  high  speed  (9600  Baud)  telephone  lines. 

The  low  speed  links  will  communicate  with  SEL  86  through  a 
program  in  the  SEL  called  TSS  (Terminal  Support  System) . Thus  the 
minicomputers  must  act  as  terminals  when  communicating  with  the 
SEL  86.  The  high  speed  links  will  communicate  with  the  SEL  86 
through  a high  speed  multiplexer  called  a MAX.  A high  speed  link 
is  a dedicated  link  and  is  active  at  all  times,  so  no  "loggin"  is 
necessary,  but  a great  deal  of  data  formulating  is  required.  The 
SEL  86  driver  for  the  high  speed  links  allows  such  features  as 
downstream  loading,  mass  data  storage  for  the  experiment  on  the 
SEL,  and  transmission  of  both  binary  and  ASCII  data. 

The  mini  and  microcomputers  in  the  laboratory  consist  of  the 
following  makes  and  models:  HP  2100  (paper  tape,  cassette,  mag 
tape,  and  disk  operating  systems) , PDP-11  paper  tape  and  floppy 
disk  system) , PDP-8  (disk  system) , Varian  620/L  (disk  system) , and 
Control  Logic  (Intel  8080)  (paper  tape  system). 
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Figure  1.  AFML  Digital  Computer  Hierarchy. 
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2.2  SELECTION  OF  A LABORATORY  STANDARD  MICROCOMPUTER 


At  the  beginning  of  this  contract  the  Materials  Laboratory 
was  using  the  Control  Logic  Microcomputer  (Intel  8080)  for  several 
ongoing  automation  experiments.  The  limits  of  this  machine  were 
becoming  very  restrictive  and  it  was  apparent  that  another  micro- 
computer should  be  chosen  for  future  projects.  The  limitations 
of  the  Control  Logic  system  included:  1)  Backplane  wiring,  2)  8 
bit  processor,  3)  limited  high  level  language  capabilities,  and 
4)  no  hardware  floating  point  arithmetic  capability.  At  the  pre- 
sent time  there  are  a great  number  of  new  microcomputers  (based 
on  minicomputers)  on  the  market.  These  machines  are  taking  ad- 
vantage of  new  Large  Scale  Integration  (LSI)  technologies  and 
offering  large  computing  power  at  very  low  cost. 

The  market  was  reviewed  and  the  Digital  Equipment  Corporation 
LSI-11  chosen  to  be  the  Materials  Laboratory  standard  microcomputer 
where  practical.  The  LSI-11  is  the  newest  in  the  PDP-11  series. 

It  is  a 16  bit  machine  with  many  DEC  options  as  well  as  non-DEC 
options.  The  LSI-11  was  chosen  because  it  has  the  following 
features: 

(1)  16  bit  registers,  memory,  I/O 

(2)  No  backplane  wiring 

(3)  FORTRAN  software  capabilities 

(4)  Floating  point  hardware  optional 

(5)  Flexible,  modular  construction 

(6)  Available  analog  interface 

(7)  DEC  software  support 

(8)  DEC  maintenance  support 

(9)  Reasonable  cost  (less  than  Control  Logic) 

At  the  time  that  the  LSI-11  was  chosen  there  were  requirements 
for  five  data  acquisition  systems,  and  since  that  time  additional 
proposals  and  preliminary  designs  using  the  I SI-11  have  been  gener- 
ated. Due  to  the  fact  that  it  is  a true  state-of-the-art  machine, 
the  LSI-11  should  remain  as  the  laboratory  standard  microcomputer 
(where  practical)  for  many  years  to  come. 
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2.3  SPECIFIC  APPLICATIONS 


One  of  the  activities  of  this  contract  has  been  the  specifica- 
tion, design,  and  implementation  of  data  acquisition  and  process 
control  systems  to  automate  the  experiments  in  the  Materials  Lab- 
oratory. This  activity  has  involved  three  computers:  Control 
Logic  (Intel  8080),  DEC  LSI-11,  and  HP  2100.  The  automation  of 
experiments  with  these  machines  is  intended  to  become  part  of  the 
Materials  Laboratory  Computer  Hierarchy  (see  Figure  1) . 

2.3.1  Liquid  Adsorption  Experiment  Automation 

The  liquid  adsorption  experiment  of  the  Surface  Inter- 
actions Branch  studies  the  rate  of  adsorption  and  desorption  of  a 
fibrous  sample  with  various  concentrations  of  solvent  and  solute 
liquids  passing  over  it.  Using  a valve  network  as  shown  in  Figure 
2 pumps  cause  various  concentrations  of  solvent  and  solute  liquids 
to  flow  into  the  sample  chamber  and  the  adsorption  of  the  sample 
is  measured  with  the  defractometer . 

The  liquid  adsorption  experiment  of  the  Surface  Anal- 
ysis Group  (AFML/MBM)  has  been  implemented  using  a Co^Jtrol  Logic 
Microcomputer.  Figures  2 and  3 show  block  diagrams  of  this  experi- 
ment and  control  system.  A Technical  Report  (AFML-TR-76-XX)  de- 
scribes the  experiment  and  automation.^  This  system  (hardware  and 
software)  was  initially  fabricated  by  AFML/DOC  personnel,  but  due 
to  changes  in  assignments  and  manpower,  the  debug,  redesign,  and 
final  implementation  (hardware  and  software)  were  done  under  this 
contract.  The  software  flowchart  is  shown  in  Figure  4.  A listing 
of  the  program  is  given  in  Appendix  A.  The  control  system  is  used 
for  two  purposes  - first  to  take  the  place  of  a human  operator  in 
overnight  experiments,  and  secondly  to  automate  the  data  gathering 
so  that  data  can  be  input  to  a data  reduction  program  in  a more 
convenient  form.  The  software  for  this  experiment  was  generated 
using  PLM,  a high  order  language  developed  by  INTEL  Corporation  for 
use  on  their  8080  microprocessor.  AFML/DOC  had  access  to  the  PLM 
compiler  on  the  GE  time-share  system.  The  program  was  written, 
stored,  and  compiled  on  the  GE  time-share  system,  and  the  cA>ject 
file  punched  on  paper  tape  to  load  into  the  Control’  Logic  System. 
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Figure  3.  Control  System  for  Surface  Analy 
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The  four  states  of  the  valves  are  called  states  A,  B, 

C,  and  the  rest  state.  These  states  relate  to  the  quantity  of 
solvent  and  solute  passed  over  the  sample.  The  rest  state  is  the 
state  in  which  the  valves  are  not  controlled  by  the  computer,  but 
rather  by  switches  on  the  front  panel.  In  state  B a solution  of 
solvent  and  solute  is  passed  over  the  sample.  In  state  C a solu- 
tion of  solvent  only  is  passed  over  the  sample.  The  solvent  to 
solute  ratio  is  stepped  in  increments  from  0 to  1 and  back  to  0. 

The  initial  and  final  states  of  the  experiment  are  state  A.  In 
state  B and  state  C the  interferometer  measures  the  adsorbtion  of 
the  sample  and  the  control  system  takes  data  from  the  interferom- 
eter, watching  for  the  signal  to  stabilize.  Once  stability  has 
been  reached,  the  control  system  changes  the  state  of  the  valves 
to  the  next  state.  Throughout  this  process  the  relative  amounts 
of  solvent  and  solute  are  changed  by  the  control  system  through 
the  stepping  motors.  In  the  past,  one  experiment  required  constant 
attention  from  an  operator  and  the  constant  taking  of  data  by  an 
operator.  Often  one  experiment  could  run  for  48-72  hours.  The 
computer  control  system  has  greatly  relieved  the  operators  task 
and  made  data  gathering  much  easier. 

2.3.2  Polymer- Surface  Interaction  Experiment  Automation 

The  polymer-surface  interaction  experiment  of  the 
Surface  Interactions  Branch  studies  the  gas  volumetric  adsorption 
properties  of  a material  sample  which  is  exposed  to  various  con- 
centrations of  gases.  The  chamber  containing  the  gases  and  sample 
is  analyzed  using  a mass  spectrometer.  This  automation  required 
the  control  of  a UTI-IOOC  Mass  Spectrometer  by  an  HP  2100  Disk 
System  minicomputer.  The  interface  was  configured  using  a Hewlett- 
Packard  universal  parallel  interface  and  a Hewlett-Packard  8 channel 
Analog  to  Digital  Converter.  The  block  diagram  of  this  automation 
system  is  shown  in  Figure  5.  The  parallel  interface  required  some 
signal  conditioning  in  order  to  drive  some  relays  in  the  mass 
spectrometer.  The  cable  and  signal  conditioner  schematics  are 
shown  in  Figure  6 . The  controls  of  the  mass  spectrometer  such  as 
operation  mode,  amplitude  range,  and  sweep  start  are  handled  by  the 
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Figure  6.  HP  2100  to  Mass  Spectrometer  Interface. 


HP  2100.  The  mass  spectrometer  has  two  analog  outputs,  one  that 
indicates  mass  number,  and  one  that  indicates  amplitude. 

The  first  step  in  the  software  development  for  this 
system  was  to  modify  the  disk  operating  system  of  the  HP  2100  to 
include  the  new  interfaces  (mass  spectrometer  and  Analog  to  Digital 
Converter) . This  also  involved  writing  an  Assembly  Language  driver 
for  the  mass  spectrometer.  This  driver  has  been  written  so  that 
it  can  be  called  from  a FORTRAN  program.  The  final  program  to 
gather  the  data  from  the  mass  spectrometer,  store  it  on  disk,  re- 
duce the  data,  and  print  it  out  is  being  written  by  AFML/DOC  per- 
sonnel . 

2.3.3  Cut-Bar  Data  Acquisition  System 

A Control  Logic  Microcomputer  was  used  to  automate  a 
Cut-Bar  experiment  of  the  Space  and  Missiles  Systems  Branch  (AFML/ 
MXS) . This  experiment  involves  a test  sample  made  of  different 
materials  to  be  heated  at  one  end  and  temperature  measurements  to 
be  taken  at  many  locations  along  the  sample  and  recorded.  The 
block  diagram  of  this  experiment  and  the  data  acquisition  system 
is  shown  in  Figure  7.  A listing  of  the  PLM  program  for  this  ex- 
periment is  given  in  Appendix  B.  Before  this  automation  it  was 
necessary  for  an  operator  to  periodically  (twice  per  hour)  manually 
record  the  voltage  on  each  thermocouple  (usually  20  in  each  experi- 
ment) to  detarm.ine  the  time  required  to  reach  steady  state  tempera- 
tures. Also,  before  the  automation  no  data  could  be  taken  at  night 
and  the  heaters  would  be  shut  off  in  the  evening  to  be  restarted 
the  next  morning.  Under  the  present  system,  the  operator  need  only 
to  initialize  the  system  (enter  the  number  of  thermocouples  and  the 
time  between  readings  on  a teletype)  and  then  the  experiment  can  be 
run  unattended. 

2.3.4  Data  Acquisition  of  Lubricants  and  Slip  Ring  Experiment 

One  of  the  projects  of  the  Liibricants  and  Tribology 
Branch  (AFML/MBT)  the  design  of  accelerated  life  tests  for  slip  ring 
bearings  that  are  used  in  satellite  applications.  Before  accelerated 


Figure  7.  Block  Diagram  of  the  Cut-Bar  Experiment. 


life  tests  can  be  designed  the  mechanisms  of  failure  must  be  under- 
stood. To  determine  the  failure  mechanisms  three  highly  instru- 
ments ted  experimental  bearing  rigs  have  been  designed.  One  rig 
will  be  used  for  long  term  studies  which  require  continuous  opera- 
tion for  times  that  may  exceed  one  year.  The  other  two  rigs  are  to 
be  used  for  relative  short  term  studies  in  which  the  slip  ring  de- 
sign, or  the  environmental  conditions  can  be  readily  changed. 

Some  of  the -parameters  that  have  to  be  measured  are  brush  speed, 
brush  force* in  three  dimensions,  pressure,  electrical  noise  gen- 
erated across  the  slip  ring,  and  environmental  constituents. 
Parameters  that  are  to  be  varied  include  brush  and  groove  geo- 
metrical design,  brush  and  groove  material  composition,  and  lubri- 
cants . 

The  data  acquisitions  systems  of  these  three  systems 
are  based  oh  the  Digital  Equipment  Corporation  LSI-11  microcom- 
puter and  are  very  similar.  The  major  differences  between  systems 
is  the  number  and  types  of  peripherals.  Thus  the  operating  sys- 
tem will  be  identical;  however,  the  peripheral  software  drivers 
will  not  be  identical  for  all  systems.  A block  diagram  of  the 
data  acquisition  system  for  the  long  term  slip  ring  experiment  is 
shown  in  Figure  8 . 

The  hardware  work  of  interfacing  the  peripherals  to 
the  LSI-11  has  been  carried  out  by  Lubricants  Branch  personnel. 
These  systems  may  contain  core  memory  for  power  failure  purposes. 
The  interfaces  that  have  been  configured  are  a Kenedy  Digital  Tape 
recorder,  a mass  spectrometer,  a multiplexer  and  controls  for  an 
Ampex  Analog  Tape  recorder,  a strip  chart  recorder,  and  various 
sensors  and  lights  within  the  experiment.  Several  signal  condi- 
tioners have  been  designed  and  build  to  modify  the  signals  of  the 
experiment  to  be  more  useful  to  the  Analog  to  Digital  Converter 
of  the  data  acquisition  system.  These  signal  conditioners  include 
filters,  peak  detectors,  integrators,  amplifiers,  and  others. 
Presently,  the  basic  hardware  of  the  first  system  has  been  install- 
ed, and  the  software  drivers  for  the  peripherals  are  being  written 
to  checkout  the  hardware.  When  these  systems  are  complete  it  will 


14 


15 


Figure  8.  Lubricants  Data  Acquisition  System  Block  Diagram. 


be  possible  for  the  data  acquisition  systems  to  take  data  in  any 
of  the  several  methods  on  a 24  hour  per  day  basis.  The  data  ac- 
quisition system  will  have  capabilities  of  taking  data  in  many 
modes,  including  time  base,  level  crossing,  and  relative  changes 
of  input  signals.  These  data  acquisition  systems  will  not  only 
replace  the  operator,  but  will  also  perform  many  duties  that  could 
not  be  done  manually. 

2.3.5  Data  Acquisiton  System  for  MTS  Materials  Testing  Rigs 

The  increased  work  load  upon  the  materials  testing  rigs 
of  the  AFML/MBC  required  a more  sophisticated  data  acquisition 
system  in  order  to  provide  a better  utilization  of  the  equipment. 

The  original  choice  by  the  personnel  of  MBC  was  a multichannel  data 
logger  with  a thermal  printer.  Reviewing  their  requirements  it  was 
determined  that  a microprocessor  based  data  acquisition  system 
could  better  satisfy  their  requirements  with  the  hardware  cost 
being  2/3  that  of  the  data  logger.  In  addition  the  increased  flexi- 
bility of  the  microprocessor  would  allow  many  more  tasks  to  be 
accomplished  that  otherwise  would  have  had  to  be  done  manually. 

The  microcomputer  data  acquisition  system  shown  in 
Figure  9 has  been  designed  and  parts  ordered  and  received.  There 
has  been  a delay  in  the  installation  of  this  system  as  the  MBC 
laboratory  is  presently  moving  and  the  new  area  is  not  yet  ready. 

The  software  for  this  system  will  involve  a time 
sharing  between  several  experiments,  multitask  operating  system. 

The  requirements  are  such  that  it  will  be  necessary  to  run  several 
distinct  experiments  at  once.  This  will  involve  starting,  stopping, 
and  changing  parameters  on  one  experiment  without  interrupting  other 
experiments.  It  will  be  necessary  to  take  data  in  several  modes: 
voltage  delta,  time  delta,  level  crossing,  and  peak  detection.  The 
software  for  this  system  will  be  written  in  FORTRAN  by  AFML/DOC 
personnel . 
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Figure  9.  Data  Acquisition  System  for  the  MTS  Testing  Rigs. 


2.3.6  Automation  of  a Stress  Intensity  Factor  Calibration 

Experiment  . 

An  LSI-11  microcomputer  was  used  to  automate  the  user 
interferometer  stress  intensity  factor  calibration  experiment  of 
AFML/LLN.  A detailed  description  of  this  experiment  was  submitted 
to  American  Society  for  Testing  and  Materials  (ASTM)  for  presenta- 
tion and  publication  at  the  Ninth  National  Symposium  on  Fracture 
2 

Mechanics.  This  experimental  technique  is  used  to  measure  crack 
growth  by  analyzing  the  diffracted  laser  light  from  a cracked  sample 
with  distinct  indentations  to  diffract  the  light.  As  the  crack 
grows,  the  interference  pattern  of  diffracted  light  passes  over  two 
photodetectors.  Using  peak  detection  software  techniques,  the 
crack  growth  can  be  accurately  measured.  A block  diagram  of  this 
experiment  is  shown  in  Figure  10.  After  peak  detection  has  occurred 
the  other  signals  are  correlated  to  the  photodetector  peaks.  For 
this  peak  detection,  it  was  found  that  floating  point  arithmetic 
would  be  necessary  to  achieve  the  desired  accuracy  and  speed.  A 
technique  for  the  peak  detection  was  derived  that  would  allow  fast 
analysis  of  the  data  without  loss  of  accuracy.  Using  this  techni- 
que, the  photodetector  signals  can  contain  frequencies  as  high  as 
2 kHz,  well  above  the  specified  requirement,  and  yet  be  analyzed 
accurately . 

The  hardware  for  this  experiment  has  been  specified,  ordered, 
recei’'cu,  and  constructed.  The  FORTRAN  software  is  approximately 
75  percent  complete  and  all  assembly  language  drivers  have  been 
written.  The  peak  detection  methods  have  been  modeled  using  FOR- 
TRAN and  found  to  be  very  accurate.  The  operating  system  for  this 
experiment  is  quite  uncomplicated  in  that  is  dedicated  to  a single 
experiment. 

2.3.7  Support  of  the  LSI-11  Software  D'  velopment  System 

The  LSI-11  software  development  system  consists  of  a 
PDP-11/03  with  12K  of  MOS  memory,  an  RXOl  dual  drive  Floppy  disk,  a 
high  speed  paper  tape  punch  (30  cps) , and  a teletype.  The  system 
runs  under  the  RT-11  operating  system  and  has  FORTRAN  IV  capabilities. 
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Figure  10.  Dlock  Diagram  of  Stress  Intensity  Factor  Experiment  Automation 


Our  usage  of  the  RT-ll/FORTRAN  system  has  been  built  from  the  four 
master  disks  provided  by  DEC.  This  includes  such  system  programs 
as  the  Editor,  the  Assembler,  the  FORTRAN  Compiler,  and  the  Linker. 
Assembly  Language  programs  (callable  by  FORTRAN  programs)  have  been 
written  to  drive  the  Analog  to  Digital  Converter,  Digital  to  Analog 
converter,  and  the  real  time  clock.  In  addition,  system  software 
has  been  written  to  drive  the  high  speed  tape  punch.  This  software 
included  modifications  to  the  RT-11  operating  system.  A listing 
of  the  PDP-11  Assembly  Language  Program  to  drive  the  high  speed 
punch  can  be  found  in  Appendix  C.  The  hardware  interface  and  cables 
were  also  configured  to  drive  the  high  speed  punch. 

2.3.8  Data  Acquisition  System  for  the  LHMEL  Facility 

The  Laser  Hardened  Materials  Evaluation  Laboratory 
(LHMEL)  facility  has  been  improved  with  the  addition  of  a dedicated 
Hewlett-Packard  9640A  Multiprogramming  System  and  a communication 
link  to  the  SEL  86  computer  system.  Software  has  been  developed 
to  acquire  experimental  data  and  to  perform  real  time  analyses  of 
these  data,  permitting  the  operator  to  rapidly  alter  experimental 
parameter  values  for  best  effect. 

The  experimental  data  that  are  acquired  are; 

(1)  Various  experimental  pre-  and  post-run  operating 
parameter  values  (entered  by  operator) 

(2)  The  system  start-up  time 

(3)  The  laser  power  (20  Hz  sample  rate) 

(4)  The  air  temperature  (20  Hz  sample  rate) 

(5)  The  static  air  pressure  (20  Hz  sample  rate) 

(6)  The  total  air  pressure  (20  Hz  sample  rate) 

(7)  The  time  the  laser  shutter  opens  and  closes 

(8)  The  time  the  burn- through  wire  in  front  of  the 
target  is  cut  by  the  laser  beam 

(9)  The  time  the  burn- through  wire  behind  the  target 
is  cut  by  the  laser  beam. 
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The  principal  component  of  the  Hewlett-Packard  system 
is  a 21MX  computer,  operating  under  the  Hewlett-Packard  supplied 
RTE-C  operating  system.  To  begin  an  experimental  run,  the  operator 
uses  RTE-C  to  activate  the  program  OZ,  and  enters  the  operating 
parameter  values  that  have  changed  since  the  last  experimental  run. 
The  operator  then  has  the  program  OZGO  activated  to  control  the 
acquisition  of  the  experimental  data.  Initialization  of  the  laser 
control  sequencer  by  the  operator  starts  the  program  data  acquisition 
sequence . 

When  the  experiment  terminates,  the  operator  uses  RTE-C 
to  activate  the  program  OXPST  to  enter  the  post-run  operating 
parameters  values.  If  desired  the  operator  can  have  the  program 
DADMP  activated  to  inspect  the  raw  data  acquired  during  the  run. 

If  the  data  are  acceptable,  the  SEL  86  computer  program  can  be 
activated.  Within  this  program  the  operator  sends  the  raw  data 
over  a communication  link  to  the  SEL  86  computer  and  initiates  the 
execution  of  the  analysis  program  SELOZ  (on  the  SEL  86) . This 
program  performs  the  following  computations: 

(1)  The  run  time  (time  of  shutter  closure  minus 
system  start-up  time) 

(2)  The  shutter  time  (time  of  shutter  closure  minus 
time  shutter  opened) 

(3)  The  burn-through  time  (back  burn-through  wire  time 
minus  front  burn-through  wire  time) 

(4)  The  laser  power  curve  over  the  run  time  interval 

(5)  The  laser  power  threshold  time  (the  time  when  the 
laser  power  becomes  less  than  (or  equal  to)  the  laser  power  thres- 
hold going  down  minus  the  time  when  the  laser  power  becomes  greater 
than  (or  equal  to)  the  laser  power  threshold  going  up)^ 

(6)  The  maximum,  minimum,  and  average  values  of  laser 
power  over  the  shutter  time  interval 

(7)  The  maximum  and  average  values  of  laser  power 
over  the  laser  power  threshold  time  interval 

(8)  The  total  energy  over  the  laser  power  threshold 

time  interval 
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(9)  The  diameter  and  area  of  the  beam  on- target 

(10)  The  laser  power  on-target  curve  over  the  run 
time  interval 

(11)  The  average  power  on- target  over  the  shutter 
time  interval 

(12)  The  average  power  density  on-target  over  the 
shutter  time  interval 

(13)  The  energy  density  on-target  over  the  shutter 
time  interval 

(14)  The  average  power  on-target  over  the  burn- 
thorugh  time  interval 

(15)  The  average  power  density  on-target  over  the 
burn-through  time  interval 

(16)  The  energy  density  on-target  over  the  shutter 


time  interval 

(17) 

interval 

(18) 

interval 

. (19) 

interval 

(20) 

(20)  The  RMS  value  of  static  air  pressure  over  the 
shutter  time  interval 

(21)  The  total  air  pressure  curve  over  the  run  time 

interval 

(22)  The  RMS  value  of  total  air  pressure  over  the 
shutter  time  interval,  and 

(23)  The  Mach  number. 


The  data  and  analytical  results  are  conveniently  for- 
matted and  printed  on  a high  speed  printer.  They  are  also  archived 
on  magnetic  tape.  CALCOMP  quality  plots  are  made  of  the  power  on- 
target,  total  temperature,  static  pressure,  and  total  pressure 
curves.  A summary  of  the  analysis  is  printed  for  the  operator  on 
the  LHMEL  teletype.  The  operator  can  use  this  summary  to  select 
the  experimental  conditions  for  the  next  run.  The  users  manual 
is  found  in  Appendix  D. 
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SECTION  III 

COMPUTER  PROGRAMS  FOR  RE-ENTRY  AND  LASER  APPLICATIONS  OF 
ABLATION/THERMAL  EFFECTS 


3 . 1 INTRODUCTION 

During  the  last  fifteen  years  many  computational  techniques 
have  been  developed  to  simulate  the  ablation/thermal  effects  occur- 
ring in  noise  tips,  heat  shields,  and  rocket  nozzles.  With  some 
modifications  these  computer  programs  can  also  be  applied  to  analyze 
the  ablation/ thermal  response  of  materials  to  laser  heating  and 
provide  guidelines  in  the  development  of  new  materials.  AFML  has 
sponsored  two  research  programs  to  select  appropriate  computational 
techniques  for  ablation/ thermal  analyses.  These  programs  have  led 
to  two  sets  of  computer  programs;  Aerotherm  Prediction.  Procedure 
for  Laser  Effects  (APPLE) , and  Unified  Ablative  Material  Thermal 
Response  Analysis  Procedure  (UNIFIED) . While  these  two  sets  of 
computer  programs  contain  many  similar  mathematical  algorithms, 
each  has  special  features  designed  for  their  particular  application. 
For  instance  APPLE  permits  in-depth  absorption  of  radiant  heating, 
while  UNIFIED  has  additional  codes  for  boundary  layer  calculations 
of  re-entry  environments  and  also  a discrete  model  capability. 

This  section  is  divided  into  two  parts;  each  part  describes 
the  application  of  one  set  of  computer  codes  to  current  AFML  prob- 
lems. The  application  of  the  UNIFIED  codes  has  been  used  to  deter- 
mine the  sensitivity  of  a one  dimensional  thermal  response  solution 
of  carbon  phenolic  in  an  arc-jet  environment  to  pertubations  of  its 
thermophysical  properties.  Specific  applications  of  the  APPLE  code 
include  the  heating  of  ZnSe  and  KCl  CO2  laser  output  coupler  win- 
dows and  simulation  of  a quartz  lampbank  heating  of  painted  alumi- 
num in  a windtunnel. 
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3.2  APPLICATIONS  OF  THE  APPLE  CODE 


3.2.1  Introduction 

APPLE  consists  of  a set  of  computer  codes  developed 

5 2 

to  describe  material  response  to  low  power  (<10  kw/cm  ) laser 
radiation  such  as  typified  by  an  infrared  CO2  laser  operating  in 
the  CW  mode.  The  basis  of  these  codes  is  the  aerospace  technol- 
ogy computer  codes  developed  for  combined  radiative  and  convective 
heating  conditions.  The  conversion  of  these  codes  to  laser  prob- 
lems is  feasible  because  low  power  heating  rate  is  similar  to  that 
experienced  by  a high  performance  re-entry  vehicle  or  a far  plane- 
tary entry  probe.  Because  these  codes  have  been  developed  speci- 
fically for  aerospace  applications  and  not  laser  heating,  a per- 
fect transition  is  not  expected  in  solving  problems  involving  laser 
heating.  Difficulties  arising  from  the  application  of  these  codes 
are  not  due.  to  numerical  computational  technique;  but  rather  to  the 
assumptions  on  thermophysical  and  thermochemical  properties, 
necessitated  by  the  lack  of  data  on  the  material  of  interest. 

The  following  descriptions  are  of  the  main  routines 

3 

comprising  APPLE  and  are  taken  from  the  APPLE  users  manual. 

The  Aerotherm  Chemical  Equilibrium  (ACE)  computer 
program  is  'an  extremely  versatile  code  for  calculating  quantities 
of  importance  to  a broad  variety  of  thermochemical  processes.  The 
thermochemical  processes  treated  may  be  divided  into  two  categories; 
closed  systems  and  open  systems.  Closed  systems  are  those  for 
which  the  relative  amounts  of  each  chemical  elements  in  the  system 
are  constant  if  prespecified.  Open  systems  are  those  for  which 
the  relative  amounts  of  chemical  elements  depend  on  various  mass 
transfer  rates  due,  for  example,  to  boundary  layer  convection  or 
solid  surface  degradation.  The  ACE  program  can  treat  both  systems 
in  chemical  equilibrium  and  systems  for  which  certain  reactions 
are  kinetically  controlled.*  The  most  important  objective  of  the 

*The  kinetic  controlled  reactions  routine  (KINET)  was  deleted  for 
use  in  the  APPLE  computer  program  library. 
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the  ACE  program  when  used  in  conjunction  with  Aerotherm's  CMA/CMAC 
and  ASTHMA  computer  codes  is  to  generate  tables  of  gas^phase  prop- 
erties at  the  wall  interface  as  a function  of  three  independent 
variables:  pressure,  dimensionless  gas  blowing  rate,  and  dimension- 
less char  blowing  rate  (p,  B',  and  B'  respectively).  In  addition 
to  calculating  the  chemical  and  thermodynamic  state  of  a variety  of 
systems,  the  ACE  program  also  calculates  and  prints  out  some  trans- 
port properties  (e.g.,  viscosity,  thermal  conductivity)  appropriate 
to  that  state. 

The  Charring  Material  Thermal  Response  and  Ablation 
Program  (CMA)  is  an  implicit,  finite-difference  computational  pro- 
cedure for  computing  the  one-dimensional  transient  transport  of 
thermal  energy  in  a three-dimensional  isotropic  material  which  can 
ablate  from  a front  surface  and  which  can  decompose  in  depth.  De- 
composition reactions  are  based  on  a three-com.ponent  model.  The 
program  permits  up  to  eight  different  back-up  materials  of  arbitrary 
thickness.  The  back  wall  of  the  composite  material  may  transfer 
energy  by  convection  and  by  radiation. 

The  ablating-surface  boundary  condition  may  take  one 
of  three  forms : 

Option  1 - General  convection-radiation  heating  with 
coupled  mass  transfer,  using  a transfer 
coefficient  approach,  including  the  effects 
of  unequal  heat-  and  mass-transfer  coeffi- 
cients (nonunity  Lewis  number)  and  unequal 
mass  diffusion  coefficients.' 

Option  2 - Specified  surface  temperature  and  surface 
recession  rate. 

Option  3 - Specified  radiation  view  factor  and  in- 
cident radiation  flux,  as  functions  of 
time,  for  a stationary  surface. 

Any  combination  of  options  may  be  used  for  a single 
computation.  Option  3 is  appropriate  to  cool  down  after  termina- 
tion of  convective  heat  input  and  is  often  useful  in  conjunction 
with  Options  1 and  2. 
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The  Axi-Symmetric  Transient  Heating  and  Material  Ab- 
lation Version  3 (ASTHMA3)  Computer  Program  is  a transient  heat 
conduction  program  for  two-dimensional,  axi-symmetric  bodies. 
Multiple  non-charring,  anisotropic  materials  may  be  studied.  The 
surface  boundary  condition  has  three  options,  including  an  unusually 
general  thermochemical  erosion  or  ablation  condition  as  well  as 
simplified  radiation  and  specified  temperature  options.  The  in- 
depth  temperature  prediction  is  of  the  familiar  explicit  finite- 
difference  type.  It  allows  a completely  general  finite  difference 
mesh  layout  relative  to  the  physical  r-z  axes , and  accounts  for 
anisotropic  heat  conduction  effects.  The  heated  surface  boundary 
condition  is  an  unusually  general  thermochemical  type.  It  accounts 
for  two  specific  kinetically  controlled  carbon  reduction  reaction, 
one  kinetically  controlled  water  gas  shift  reaction,  and  any  num- 
ber of  gas-phase  equilibrium  reactions  for  any  combination  of  ab- 
lating materials  and  environments.  Chief  applications  for  the 
computer  program  are  rocket  nozzles  and  entry  vehicles.  The  follow- 
ing sketch  serves  to  clarify  various  physical  aspects  of  the  abla- 
tion problem  treated. 


Computation 

Calculate  radiation  flux  to  wall 

-Determination  of  boundary  layer 
edge  state 

-Account  for  boundary  layer  trans  - 
port  of  mass  and  energy  to  and  from 
surface 

-Determine  chemical  state  at  surface 
and  perform  energy  balance 
coupling  to  in-depth  solution 

““In -depth  thermal  response 
calculation 


Physics 


boundary 
layer  edge 


mam 
stream 


k. 
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The  JANAF  routine  is  a special  purpose  code  that  has 
been  incorporated  into  the  COUPLE  computer  program.  It  assembles 
a molecular  species  thermochemical  data  (JANAF)  file  in  the  appro- 
priate format  for  input  into  the  ACE  program.  JANAF  performs  these 
functions  by  searching  in  the  JANAF  data  tape  for  the  molecular 
species  that  can  be  formed  from  a prescribed  list  of  atomic  elements 
that  is  specified  in  the  input.  . 

The  COUPLE  preprocessor  routine  is  a user  oriented 
program  specifically  developed  to  execute  singly  or  coupled  a 
variety  of  Aerotherm  surface  state  thermochemistry  and  in-depth 
ablation  and  thermal  response  computer  programs.  It  generates 
appropriate  input  data  for  the  computer  programs  that  are  to  be 
executed  from  simple  execution  directives  and  very  little  input 
data.  Version  1 couples  the  ACE  and  JANAF  routines  to  obtain  sur- 
face state  thermochemistry  solutions  that  are  used  as  lan  input  to 
the  CMA  and  ASTHMA  computer  programs.  Version  2 couples  the  ACE, 
JANAF,  and  CMA  routines,  while  version  3 couples  ACE,  JANAF,  CMA, 
and  ASTHMA  routines.  j 

While  running  APPLE  coupled  decreases  total  turn  around 
time  for  analysis  as  only  a single  job  is  submitted,  it  is  usually 
desirable  to  run  APPLE  uncoupled.  Thus  no  parameter  values  are 
assigned  by  default  and  the  validity  of  assumptions  made  can  be 
more  readily  determined.  In  addition,  if  the  problem  doesn't  run 
correctly,  debug  is  easier  as  the  responsible  routines  can  more 
easily  be  located.  Brief  descriptions  of  problems  solved  with  the 
APPLE  code  under  this  contract  are  included  below.  A, more  detailed 
description  of  these  problems  can  be  found  in  a series  of  AFML  tech- 
nical memos^  and  in  a technical  report  by  Rondeau  and  Ford.^^ 

3.2.2  Thermal  Response  of  "Painted"  Aluminum  to  Radiant 

Heating 

Currently,  various  elastomeric  coatings  on  aluminum 
plates  mounted  in  a windtunnel  are  being  exposed  to  radiant  heating 
by  a quartz  lampbank.  The  objective  of  this  experiment  is  the  as- 
sessment of  the  nuclear  flash  protection  afforded  by  the  elasto- 
meric coating  applied  to  an  aircraft's  outer  skin.  A model  of 
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of  this  experiment  consisting  of  the  thermophysical  properties  of 
a polyurethane  coating  material  and  aluminum  2024T3,  the  physical 
dimension,  and  the  boundary  conditions  (heating  rate  and  heat  trans- 
fer coefficient)  were  used  with  the  APPLE  code  to  calculate  the  back 
face  temperature  of  the  aluminum  substrate  as  a function  of  time. 

The  agreement  with  experimentally  obtained  back  face  temperature 
5 

was  quite  good.  The  use  of  APPLE  does  not  supplant  the  need  to 
do  the  experiment  but  rather  compliments  it.  The  experiment  helps 
determine  physical  model  to  be  used  as  input  to  APPLE,  i.e.,  knowl- 
edge of  whether  the  coating  would  peel,  ablate,  or  blister  during 
the  experiment.  Once  agreement  between  experimental  and  calculated 
results  is  obtained,  the  APPLE  code  can  then  be  used  to  predict  the 
results  for  extrapolated  experimental  conditions. 

In  addition  APPLE  was  used  to  parametrically  study  the 
effect  of  coating  thermal  conductivity  on  the  surface  and  in-depth 
temperatures  of  a polyurethane  coated  aluminum  sample  subject  to  a 
thermal  pulse  in  a wind tunnel.  The  results  show  that  temperature 
changes  from  a +20%  variation  in  the  thermal  conductivity  of  the 
coating  are  relatively  minor. 

The  code  was  also  used  to  model  an  in-flight  scenario 
where  the  wing  section  of  an  aircraft,  traveling  at  Mach  0.85  and 
30,000  feet  is  exposed  to  the  same  heat  flux  as  that  produced  in 
the  laboratory  experiment.  Comparisons  between  the  flight  simu- 
latif  n ar.d  laboratory  experiments  have  been  presented  along  with 
the  parametric  study  of  the  coating  thermal  conductivity  in  Refer- 
auL-es  V cind  9. 

3 , ? 3 ' Therma 1 and  Stress  Analysis  of  the  AFWL  Zinc  Selenide 

La^er  VTinj^w 

Operating  conditions  of  the  Air  Force  Weapons  Labora- 
" V n 'trie  Di3c.heT.gr;  Coaxial  Laser  (EDCL)^^  were  used  to  calcu- 
1 rs  ;»>  H eynmit'tric  (21'-)  temperature  profiles  in  the  zinc  sele- 
nide output  coupler  CO^  laser  windows.  These  temperatures  were 

uaod  to  oalouioto  Ik  , diexmally  induced  otresses  in  the  win- 
dow vre  < f ; a .Tf.i.'  ..<■  • j.  c-  wo^l  w ithin  design  expectations 

for  th  ; ■:  ' . •> 


3.2.4 


Computed  Temperature  Response  cf  the  AFML  Flat-Top 
Laser  Window 

AFML's  10  Kilowatt  Flat-Top  laser  is  based  on  the  design  of 
AFWL  EDCL  laser. Because  of  the  greater  absorptance  of  the  AFML 
zinc  selenide  windows,  as  determined  by  post-mortem  examination, 
sufficiently  high  thermally  induced  stresses  resulted  in  two  win- 
dow failures  under  normal  operating  conditions.  The  APPLE  codes 
were  used  to  verify  that  thermally  induced  stresses  were  a possible 
failure  mechanism  and  to  generate  curves  of  maximum  thermally  in- 
duced stress  versus  time  for  various  power  levels.  These  calcula- 
tions provide  safe  operating  criteria  to  prevent  failure  of  the 
laser  coupler  window.  This  work  is  described  in  detail  in  Re- 
ference 8. 

3.3  UNIFIED  COMPUTER  PROGRAMS 

3.3.1  Introduction 

The  primary  objective  in  Unified  set  of  computer  pro- 
grams as  developed  by  Aerotherm  is  to  apply  the  advanced  computer 
technology  developed  over  the  past  15  years  (for  the  ablation  de- 
sign analysis  problem)  to  the  far  more  difficult  problem  area  of 
materials  development  guidance.  The  relative  payoffs  of  altering 
the  chemical  composition  of  a composite  material,  reductions  of 
density  and  thermal  conductivity,  and  changing  the  composites  con- 
figuration both  on  a macroscale  and  a microscale  are  all  important 
areas  that  can  be  addressed  with  the  Unified  set  of  codes.  The  one 
dimensional  computer  program,  CMA,  primary  use  has  been  to  deter- 
mine the  sensitivity  of  material  response  in  an  arc  jet  to  changes 
in  the  thermophysical  properties  of  carbon  phenolic.  This  data 
coupled  with  the  relative  uncertainties  of  the  individual,  thermo- 
physical properties  illustrate  where  the  most  beneficial  use  of 
scarce  resources  can  be  used  to  reduce  the  uncertainty  of  the 
parameters  that  most  effect  the  material  response.  Similar  analy- 
ses can  be  done  for  a re-entry  environment  and  using  the  two  dimen- 
sional computer  program  ASTHMA4;  however  the  more  interesting 
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results  will  be  using  ASTHMA4  to  obtain  material  responses  due  to 
material  configuration.  The  first  study  using  CMA  is  described 
below. 

3.3.2  Sensitivity  Study  Using  CMA 

A Carbon  Phenolic  button,  backed  by  a stainless  steel 

plate,  placed  in  an  aic  jet  was  modeled  for  analysis  with  the  CMA 

computer  code.  Fourteen  nodes  were  placed  in  the  Carbon  Phenolic 

button  and  two  nodes  in  the  stainless  steel  backup  plate.  The 

environraenbal  conditions  used  were  those  similar  to  a typical  arc- 

jet  experiniental  test  run  (i.e.,  pressure  = 1.07  atm,  heat  input 
2 

of  1050  BTU/(ft  -sec)  for  30  seconds  and  0 for  the  remaining  70 
seconds) . The  baseline  thermophysical  properties  of  carbon  phenolic 
that  were  used  are  shown  in  Table  I.  The  properties  that  were 
varied  +10%  each  are  shown  in  Table  II. 

The  response,  as  a function  of  time,  of  the  carbon 
phenolic  material  was  measured  by  seven  parameters.  These  are  the 
interface  t;eniperatures  between  the  carbon  phenolic  and  the  stain- 
less steel  and  6 dimensionless  variables  shown  in  Table  III. 


Thirty-five  computer  runs  were  made  (2  runs  as  each 
parameter  was  varied  +10%  and  the  baseline  case) . Since  the 
material  response  functions  as  defined  in  Table  II  are  not  calcu- 
lated directly  by  the  CMA  computer  program,  the  CMA  computer  pro- 
gram was  modified  to  create  a new  file  which  contained  all  quan- 
titxef  necessaij  to  calculate  these  response  functions.  This  file 
f f.<  r.alogusd  and  then  at  a later  time  used  as  an  input  file  to 
H jter  prcgraiii  RES  which  then  calculated  the  individual  re- 

sponse lunction.  Ihis  computer  program  is  listed  in  Appendix  E. 


The  -ji-nsitivi  by  coefficients,  defined  as 


AR. 


AX. 


where 


t)  J 

/■  R -K  *-,hr  -•■nf>nrre  In  the  j material  response  function  due  to 

a /X.  charge  made  to  the  material  input  parameter,  are  deter- 

n.nel  oy  the  computer  codes.  The  sensitivity  coefficients  for  R^ 

'^he  interface  tempera+ere)  and  the  experimental  uncertainty  in 

12 

ttie  matin'i-jl  s i > t-  .•  ' I:'."' ted  in  Table  TV.  The 
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TABLE  I 

thermophysical  properties  of  carbon  PHENOIIC 
(CCAZ  CLOTH/5C  1008  RESIN) 


PYROLYSIS  GAS  ENTHALPY 


-1899.0 

-761.9 

8868.6 

2894.7 

3953.5 

5747.2 


INTERNAL  DECOMPOSITION  KINETIC  DATA 

i 

"vi 

B) 

Ej/°R 

Ib/ft^ 

Ib/ft^ 

sec'^ 

°R 

1 (resin) 

19.34 

0. 

1.4x10^ 

1.54x10^ 

2 (resin) 

58.00 

29.24 

4.8x10’ 

3.68x10^ 

3 (reinforcement) 

98.55 

98.55 

0. 

0. 
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MATERIAL  RESPONSE  PARAMETERS  (R|) 


L 


dx 


0 

r^dx 

•'o 


M • 


fp  ^ 

f ' ^in 

» r\ 


Tgp  - INTERFACE  TEMPERATURE 
k - THERMAL  CONDUCTIVITY 
P~  DENSITY  OF  THE  iVlATERIAL 
Cp  - SPECIFIC  HEAT  OF  THE  MATERIAL 
X - THICKNESS  OF  THE  MATERIAL 
qjp-  RATE  OF  HEAT  ABSORBED  BY  THE  MATERIAL 
rate  of  heat  LOST  BY  THE  MATERIAL 
hg  - ENTHALPY  OF  THE  GASEOUS  MATERIAL 
riig  - RATE  OF  MASS  LOSS  OF  THE  MATERIAL 
- ENTHALPY  OF  THE  SOLID  MATERIAL 
H^  - HEAT  OF  FORMATION 
t - TIME. 
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TABLE  IT 

SIMS!T1VIT\  COEFFICIENTS  FOR  INTERFACE  TEMPERATURE 


PYROLYSIS  REACTION  CONSTANTS,  RESIN  1 0.1275 

PYROLYS I S REACTION  CONSTANTS,  RES  IN  2 0. 2270 

HEATOFFORAAATION  0.1287 


total  uncertainty  of  a material  response  function  AR^is  given  by 


AR.  = 
D 


AR. 


AX. 

1 


AX. 

1 


(1) 


AR. 


with  -r^  bieng  the  sensitivity  coefficient  and  AX.  the  experimental 
uncertainty  of  the  i material  input  parameter.  An  appreciable 
reduction  in  AR^  can  be  made  by  reducing  the  uncertainty  AX^  only 
if  the  product  involving  in  Equation  (1)  is  large. in  comparison 
to  other  terms  in  the  sum.  The  parameters  that  most  significantly 
effect  the  uncertainty  of  the  interface  temperature  as  determined 
by  this  study  and  listed  in  Table  IV  are: 


(1) 

CHAR  CONDUCTIVITY 

- MOD 

SENS/HIGH 

UNC 

(2) 

VIRGIN  CONDUCTIVITY 

- MOD 

SENS/HIGH 

UNC 

(3) 

CHAR  FABRIC  DENSITY 

- HIGH 

[ SENS/MOD 

UNC 

(4) 

VIRGIN  FABRIC  DENSITY 

- HIGH 

[ SENS/MOD 

UNC 

(5) 

VIRGIN  SPECIFIC  HEAT 

- MOD 

SENS/MOD  UNC 

(6) 

REACTION  CONSTANTS 

- MOD 

SENS/MOD  UNC 

Similar  analyses  are  being  conducted  for  the  other 
response  characteristics  and  are  to  be  published. 
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SECTION  IV 

COMPUTER  ANALYSIS  OF  MICROSTRUCTUPES 

4.1  COMPUTER  RECOGNITION  OF  TITANIUM  ALLOW  MICROSTRUCTURES 
‘1.1.1  Ir.trodurtion 

The  objective  of  this  project  is  to  predict  the  mechan- 
ical properties  of  titanium  alloys  by  computer  examination  of  photo- 
micrographs of  their  microstructures.  Prior  to  this  contract,  cora- 

13 

puter  codes  were  developed  to  effect  the  solution.  The  results, 
however,  were  inconclusive;  the  correlations  obtained  did  not  con- 
sistently provide  accurate  predictions  of  the  mechanical  properties 
of  samples  processed.  The  previously  developed  solution  was  basi- 
cally three  steps. 

1.  A digital  image  of  a microstructure  is  produced 

by  digitizing  a photomicrograph  negative.  The  digital  image  consists 
of  light  intensity  readings  (pixels)  taken  at  40  micron  intervals 
of  the  negative.  The  pixel  (picture  points)  intensities  are  valued 
0-255. 

2.  The  digital  image  is  operated  on  by  computer  soft- 
ware, the  particles  of  the  microstructure  are  extracted  and  their 
aeometric  characteristics  determined. 

3.  These  geometric  characteristics  of  each  image  are 

. jr -e 1 ? ' e i .with  the  measured  mechanical  properties  of  the  respective 
samples,  yielding  a relation  that  will  predict  mechanical 
propertres  from  the  geometric  characteristics  of  the  particles  in 
the  mirtoat;  iicture  of  nn  alloy  sample.  Work  performed  on  this  pro- 
jfC*  US'.’  ■ r'>''tract  was  concerned  primarily  with  improvements 

■-  LjTip’.'.ter  prccessxng  involved  in  Part  2 above. 

i . i . z Inxciai  Work 

The  step  taken  in  working  on  the  project  was  to 

become  familiar  with  the  problem  and  gain  a working  knowledge  of 
j-i.f-  rf  ' h-  ■ t < pvoc^ss  the  digital  images. 


Soon  after  work  was  underway  on  the  project,  a new 
set  of  photomicrographs  of  the  Ti-alloy  surfaces,  which  were  of 
higher  quality  than  photomicrographs  used  in  the  previous  work, 
became  available.  It  was  theorized  that  processing  these  would 
yield  better  correlation  results,  and  the  processing  of  the  new 
images  was  begun.  Several  of  the  programs  that  performed  initial 
image  processing  were  found  to  be  inoperative.  It  was  determined 
that  the  majority  of  the  problems  encountered  in  the  attempt  to 
execute  these  programs  were  due  to  changes  in  the  operating  system 
of  the  CDC  6600  computer  that  had  been  made  since  the  programs 
were  written  and  last  used  in  a production  mode.  The  computer 
codes  were  modified  to  make  them  compatible  with  the  pre.ant  CDC 
6600  operating  system  and  the  processing  of  the  new  images  con- 
tinued. 

4.1.3  Improvements  of  the  Digital  Processing 

A close  examination  of  the  available  intermediate  re- 
sults of  the  processing  done  during  previous  work  on  the  project 
was  undertaken  in  an  effort  to  gain  insight  into  possible  improve- 
ments that  could  be  made.  Conclusions  reached  were  that  image  thres- 
holding and  particle  separation/extraction  should  be  the  areas  of 
refinement.  The  digital  image  thresholding,  in  which  ‘the  raw  image 
pixel  values  of  0-255  are  reduced  to  a value  of  0 or  1,  is  the 
most  critical  stage  of  the  processing;  improvements  in  this  process 
would  cause  significant  improvement  in  the  accuracy  of  the  particle 
geometric  characteristics  determined  at  the  final  stage  of  the  image 
processing.  In  the  particle  separation/extraction  phase,  it  was 
observed  that  frequently  particles  that  were  very  close  together 
in  the  photomicrograph  were  extracted  from  the  digital  image  and 
treated  as  one  large  particle.  This  was  due  to  higher  density 
levels  of  the  pixels  in  the  area  of  the  junction  of  these  particles, 
caused  by  light  scattering  and  the  resultant  noise  when  the  photo- 
micrograph was  digitized. 
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The  two  operations  of  image  thresholding  and  particle 
separation,  however,  are  not  totally  independent  of  each  other.  It 
was  seen  that  the  accuracy  of  particle  separation/extraction  would 
be  greatly  facilitated  by  a more  accurate  thresholding  of  the  digital 
images.  The  approach  taken  to  the  improvement  of  the  image  process- 
ing was  the  development  of  a more  accurate  thresholding  process, 
and  the  later  development  of  a more  precise  method  of  particle  sepa- 
ration/extraction, if  necessary. 

4.1.4  Improvement  of  linage  Thresholding 

The  initial  step  in  establishing  an  improved  threshold- 
ing process  was  the  collection  of  data  of  raw  image  characteristics. 
This  data  consisted  primarily  of  computer  printouts  of  the  digital 
images  and  corresponding  plots  of  pixel  density  versus  frequency 
for  a selection  of  photomicrographs. 

In  the  previous  method,  the  threshold  density  D,^ 
was  determified  by 


^T  = 

where  D was  the  density  below  which  five  percent  of  the  pixels  of 

La 

the  image  were  contained,  and  D„  the  density  above  which  five  per- 

n 

con^  of  the  pixels  were  contained.  (This  represents  the  standard 
technique  for  thresholding  used  in  image  analysis.)  Through  com- 
t-ur  ;.sori  ^the  oiiginai  photomicrograph,  the  printout  of  the  raw 
' 5 tiK  <.'ensit/  versus  frequency  plot,  and  the  printout  of  the 
-’..1  uf ’".oldod  image  “or  ea^'h  sample  used,  it  was  concluded  that  this 
”■>■£1100  dj.d  uot  ccr  yield  the  best  thresholded  image. 

Several  al  tp’'-iate  thresholding  methods  were  examined. 
The  reot/o ques  evajuated  first  were  differential,  line-by-line 
thi eeaolding  methods,  and  failed  to  produce  accurately  thresholded 
The  po-ja.  ^eeuits  obtained  are  believed  to  be  due  to  over- 
all density  variation  across  the  photomicrograph,  relatively  high- 
denr-ity  a '-i  / ' themselves,  and  general  in- 

accuracies a£,?ocLotfc'>l  > jtr<  ;;  !.•>  lao  ore  dimensional  operations  to 


The  thresholding  method  arrived  at  was  essentially  a 
modification  to  the  original  method.  The  thresholding  density 
was  determined  as 

""t  = *°I 

where  and  remained  unchanged  from  the  original  method,  and 
Dj.  is  the  percentage  of  the  remaining  image  pixels, to  be  consider- 
ed as  being  within  a particle.  D^.  is  an  input  parameter  to  the 
thresholding  program.  Using  a carefully  selected  value  of  D^., 
determined  by  examining  the  density  histogram  for  each  image  proc- 
essed, the  respective  thresholded  images  were  improved,  significantly. 

4.1.5  Improvements  in  Particle  Separation/Extraction 

During  the  processing  of  the  new  photomicrographs,  the 
particle  separation/extraction  portion  of  the  computer  program  was 
found  to  be  totally  inoperative.  Due  to  the  program's  complexity 
and  lack  of  documentation,  it  was  necessary  to  examine  the  entire 
code  in  detail  in  an  effort  to  locate  and  correct  the  error (s) . 

In  doing  so,  it  was  discovered  that  numerous  major  programming 
errors  were  present;  that  in  fact  when  operative  the  codes  would 
not  perform  the  operations  as  described  in  the  furnished  docxomenta- 
tion.  Because  of  these  conditions,  the  design  and  development  of 
new  particle  separation/extraction  software  was  required.  Also, 
the  magnitude  of  the  errors  was  certain  to  have  adversely  affected 
the  results  previously  obtained  using  these  codes. 

4.1.6  Development  of  New  Particle  Separation/Extraction 

Software 

Several  methods  of  particle  separation/extraction 
were  investigated,  in  order  to  select  the  most  accurate  one  for 
implementation.  The  first  methods  investigated  employed  circular 
and  rectangular  "windows"  of  various  sizes  which  were  "moved"  over 
the  raw  digital  images;  while  examining  the  portions  of  the  image 
contained  in  the  window.  Particle  edges  are  detected  using  this 
method  by  looking  for  large  density  variations  occurring  within 


the  window.  Once  such  a density  change  is  found,  the  window  is 
then  moved  slowly  around  the  particle  by  maintaining  the  density 
differential  interior  to  the  window  while  shifting  its  position  on 
tiie  image,  thus  outlining  a particle.  This  approach  failed,  however, 
due  to  the  wide  range  of  particle  sizes,  making  the  selection  of 
the  proper  window  ar.?a  by  the  computer  program  extremely  difficult 
at  best.  (Often  a window  would  entirely  enclose  a small  particle 
in  one  part  of  an  image,  and  then  be  entirely  enclosed  itself  by 
a large  particle  in  a nearby  option  of  the  image.)  Separation  of 
two  seemingly  connected  particles  was  unpredictable  using  this 
method,  due  to  the  varying  sizes  of  these  "bridges"  connecting  the 
particles  in  the  digital  image. 

The  method  of  particle  separation/extraction  chosen 
to  be  implemented  employs  the  same  basic  concepts  as  the  original 
theory. In  addition,  the  software  incorporates  the  following: 

1.  input  parameters  should  provide  the  capability  of 
selecting  the  exact  separation  criteria  for  the  case  of  erroneously 
connected  particles; 

2.  actual  separation  should  not  be  biased  with  re- 
spect to  particle  location  in  the  image  (previous  method  was  biased 
in  an  upper-left  sense,  due  to  processing  the  image  left  to  right 
and  top  to  bottom) ; 

3.  the  software  should  be  much  more  efficient  in 
»-orTr/!  of  central  memory  required  than  its  predecessor; 

4.  intermediate  states  of  the  image  during  separation/ 

I ; .:a  .'  ill  should  De  available  for  precise  evaluation  of  software 

' O..I'  . 

The  sepa ratio  V'^r^tac^- ion  operation  was  broken  into  two  main  steps, 

Ic  p.if'.riiea  two  individual  programs: 

1.  Generation  of  interior  particle  pixel  values  (1,2, 

i.  . jpoi'  tue  location  of  a pixel  with  respect  to  the 

outer  boundary  of  the  particle) . 

2.  por'-'r- r Td  extraction  of  the  particles. 


4(' 
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The  design,  development,  and  testing  of  a program 
(PARSE?)  to  perform  the  first  step  has  been  completed.  This  program 
requires  much  less  central  memory  for  its  execution  (55,000g  words 
versus  165,000g  words  for  the  previous  program) . This  improvement 
was  made  possible  by  maintaining  a much  smaller  portion  of  the  image 
in  central  memory  during  processing.  In  addition,  a utility  progreuti 
(IM6PRN)  was  written  to  print  the  image  during  and  after  the  proc- 
essing for  software  evaluation  purposes. 

The  programming  of  the  second  step  to  process  the  image 
output  by  PARSE?  (separation  and  extraction)  was  divided  into  five 
parts:  i 

1.  separation  and  extraction  of  particle. cores  (3- 
valued  pixels) ; 

I 

2.  determination,  separation  (when  necessary),  and 
extraction  of  the  upper  particle  extremes  {2-valued  pixels) ; 

3.  determination,  separation  (when  necessary),  and 
extraction  of  the  lower  particle  extremes  (2-valued  pixels) ; 

4.  determination,  separation  (when  necessary),  and 
extraction  of  the  remaining  left  and  right  particle  extremes  (2- 
valued  pixels) ; and 

5.  reconstruction  and  output  of  the  particles. 

The  determination  and  separation  operations  stated 
above  are  effected  by  moving  a window  of  dimensions  1 pixel  by  1 
pixel  around  in  the  stated  area  of  a particle.  This  method  locates 
particles  from  their  interiors  radially  outward.  In  the  case  of 
seemingly  "connected"  particles,  the  user  may  specify  the  exact 
location  at  which  these  particles  will  be  separated.  »This  separa- 
tion actually  talces  place  at  the  user-specified  relative  location 
between  the  particle  cores  (3-valued  pixels).  For  example,  the 
user  may  specify  that  "connected"  particles  be  separated  at  the 
location  halfway  between  their  cores.  Software  to  perform  parts  1, 
2,  and  3 has  been  designed,  developed,  and  tested.  Software  to 
perform  steps  4 and  5 has  been  designed  and  is  presently  in  the 
development  stage. 


1 
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4.1.7  Additional  Software  Developed  for  Testing  and  Evalua- 
tion Purposes 

In  the  process  of  testing  and  evaluating  both  previous- 
ly existing  ana  newly  developed  software,  the  following  programs 
were  written: 

1.  PXCTlfRE  - A program  to  print  a "picture"  of  the 
raw  digital  images  using  32  grey  levels,  obtained  by  overprinting 
of  characters,  for  the  purpose  of  quick,  easy  comparison  of  the 
raw  digital  image  to  the  photomicrograph  used  to  create  it. 

2.  PARPLOT  - A program  to  reconstruct  and  graphically 
display  a microstructure  from  the  geometric  characteristics  output 
by  the  MEAD  program  PROPERT. 

3.  SAMPLOT  - A program  to  produce  plots  of  frequency 
versus  the  log  of  area  for  particles  extracted  from  images  in  pre- 
vious work  On  the  project. 

4.  FREQPLT  - A program  to  produce  plots  of  the  image 
density  histograms  produced  by  the  MEAD  program  HISTCDC. 

In  addition  to  the  work  described  above,  other  work 
perforiTied  in  association  with  this  project  consisted  of  the  follow- 
ing ; 

1.  Preparation  of  10  image  tapes  used  by  the  PAR  Cor- 

14 

poration  in  work  related  to  this  project. 

2.  Familiarization  with  the  use  of  the  OLPARS  (On- 

I i Pattern  Analvsis  and  Recognition  System)  of  the  Rome  Air  De- 
velopment Center  at  Griffiss  Air  Force  Base  for  the  purpose  of 
determining  geometric  characteristics-mechanical  properties  corre- 
l.'tio;'*!  '’’his  system  was  determined  suitable  to  perform  correla- 
c na'.ys  £ wnen  the  new  image  geometric  data  becomes  available. 


4.2  COMPUTER  MODELING  OF  CRACK  PROPAGATION  IN  THE  SURFACE  TI-ALLOYS 


4.2.1  Problem  Description 

In  work  performed  with  AFML/LLS,  an  effort  to  develop 
computer  codes  to  model  crack  propagation  across  the  surface  of 
Titanium  alloys  was  initiated  as  part  of  a continuing  effort  in  this 
area.  Given  an  alloy  surface,  the  project  objective  is  to,  use  the 
computer  to  predict  the  nature  of  a crack  growing  across  it.  At  the 
present,  this  alloy  surface  is  presented  as  an  enlarged  photograph 
of  an  alloy  surface  in  which  the  features  affecting  crack  propaga- 
tion (grain  structure,  large  particles,  etc.)  are  easily  discern- 
able  to  the  naked  eye.  Crack  propagation  is  to  be  effected  accord- 
ing to  propagation-surface  conditions  relations  specified  by  the 
software  user. 


4.2.2  System  Requirements 

Initial  analysis  of  the  problem  yielded  the  following 
requirements  for  the  developed  software  for  the  modeling  process. 


1.  Capability  of  accurate  representation  of  the  alloy 

surface . 

2.  Capability  of  easy  modification  to  that  surface 
representation . 

3.  A high  degree  of  flexibility  in  allowing  for  the 
specification  of  crack  behavior  criteria  with  respect-  to  surface 
conditions,  yet  perform  in  strict  accordance  to  that  criteria. 


Due  to  the  desirability  of  a high  degree  of  user- 
software  interaction,  the  software  developed  functions  in  an  inter- 
active mode.  There  existed  the  need  for  graphics  capability,  for 
both  establishing  and  modifying  the  surface  representation,  and  for 
displaying  intermediate  and  final  results  of  the  modeling  process. 
AFML's  TEKTRONIX  4014  graphics  terminal  and  4954  graphics  tablet 
in  conjunction  with  the  INTERCOM  system  of  the  CDC  6600  was  select- 
ed to  be  used  since  this  configuration  satisfies  the  interactive 
graphics  capabilities  required. 


4.2.3  Computer  Representation  of  the  Alloy  Surface 

The  photograph  of  an  alloy  surface  to  be  processed  is 
used  to  establish  the  digital  surface  representation  in  the  com- 
pute: . This  consxsts,  basically,  of  five  steps: 

1.  Prior  to  execution  of  the  modeling  program,  the 
boundaries  of  areas  of  homogeneous  grain  orientation  on  the  photo- 
graph are  manually  outlined.  Within  each  outlined  area,  a line 
determining  the  grain  orientation  to  be  used  for  that  area  is  drawn. 

2.  Each  boundary  and  respective  grain  orientation  line 
on  the  photograph  are  digitized  using  the  TEKTRONIX  4954  graphics 
tablet. 

3.  The  output  from  the  digitizer  is  stored  in  the 
computer.  The  computer  surface  representation  then,  in  essence, 
consists  of  a set  of  lines. 

4.  The  surface  representation  is  displayed  graphically 
on  the  terminal  screen.  Any  bad  points  resulting  for  erroneous 
digitized  input  are  corrected  by  editing. 

5.  The  final  step  is  to  eliminate  gaps  between  points 

on  the  boundary  lines,  to  insure  that  a crack  cannot  pass  between 
two  points  on  a boundary  line  and  not  be  recognized  by  the  computer 
as  ha\ing  intersected  that  boundary  line.  This  is  done  by  convert- 
ing boundary  line  point  values  to  integers  and  interpolating  points 
betw'f n the  original  points  comprising  a boundary  line  such  that 
for  consecutive  points  on  a boundary  line  (X^^,  , (X^,  Y2)  » 


ir-J 


iX^  - Xj^j  = 1 or  0, 


1^2  ' '^1'  = ^ 0* 


tr'  anv  one  ( f the  above  relations  must  have  a value  of  1. 


Tht  at'o ’c  orrrations,  with  the  exception  of  1,  are 
performed  by  programs  of  the  CRAKMOD  software  system. 


4.2.4  Linkage  Table  Generation 


Once  the  user  is  satisfied  with  the  established  sur- 
face representation,  the  table  of  point-to-point  links  used  to 
determine  crack  propagation  is  generated  by  CRAKMOD.  At  the  present 
stage  of  development,  these  links  are  determined  solely  from  the 
grain  orientations  of  each  boundary.  (The  linkage  table  structure, 
however,  is  designed  to  provide  for  additional  links  to  be  deter- 
mined by  additional  crack  propagation-surface  condition  criteria, 
as  these  relations  become  known.) 

4.2.5  Crack  Propagation 

The  propagation  by  the  computer  of  a crack  across  the 
represented  surface  is  initialized  by  the  CRAKMOD  user  specifying 
its  starting  location.  The  software  then  determines  the  path  of 
the  crack  by  chaining  together  the  points  specified  by  the  linkage 
table . 

The  main  output  from  CRAKMOD  is  the  graphical  display 
of  the  surface  and  the  crack  grown  across  it.  (Optional  display 
capabilities  include  magnified  views  of  one  or  more  user  selected 
surface  areas,  with  or  without  the  segments  of  the  crafck  (if  any) 
across  them  displayed.) 

At  the  present  stage  of  development,  the  CRAKMOD  user 
may  effect  different  crack  configurations  by  specifying  changes  in 
the  starting  location  of  the  crack  and  the  grain  orientations  for 
the  areas  enclosed  by  the  boundary  lines.  As  other  surface  condi- 
tion/crack propagation  relations  are  recognized,  all  that  will  be 
necessary  to  modify  CRAKMOD  for  their  implementation  will  be  the 
addition  of  program  modules  to  generate  the  proper  links,  and 
provide  for  the  input  of  any  additional  necessary  features  from 
the  photograph  (via  digitizer  or  terminal) . CRAKMOD  is  completely 
modular  in  structure,  such  that  the  implementation  of  the  above 
features  will  require  a minimum  amount  of  modification  to  the  ex- 
isting software. 


4.2.6  Additional  Capabilities 


In  addition  to  the  generation  of  cracks  across  the  sur- 
face representation,  CRAKMOD  is  capable  of  processing  a crack  digi- 
tized from  the  alloy  surface  photograph.  Thus,  crack  generated  by 
CRAKMOD  may  be  compared  to  an  actual  crack  in  the  surface  processed. 

There  is  the  capability  of  measuring  the  length  of  the 
crack  very  accurately;  much  more  accurately  than  is  possible  by 
measuring  directly  from  a photograph  or  an  actual  sample,  because 
the  measurement  may  be  made  along  the  length  of  the  actual  crack 
using  CRAKMOD,  rather  than  along  a straight  line  from  its  initial 
point  to  its  end. 

A listing  of  all  FORTRAN  programs  developed  for  com- 
puter analysis  of  microstructures  is  enclosed  in  Appendix  F. 


SECTION  V 

GENERAL  USER  PROGRAMS  AND  EDUCATION 

1 

5 . 1 INTRODUCTION 

A number  of  general  user  programs  were  written  for  use  by  AFML 
personnel.  These  include  programs  designed  for  graphics,  digitizing, 
data  analysis,  support  of  the  LSI-11  microprocessors  and  support  of 
the  Control  Logic  microprocessors.  In  addition,  classes  for  AFML 
CDC  6600  and  LSI-11  users  were  held  to  improve  the  utilization  of 
these  computers. 

5.1.1  General  User  Programs  for  the  CDC  6600 

DIGIPLT  is  an  easy-to-use,  completely  interactive 
software  package  that  provides  the  capability  of  digitizing,  dis- 
playing in  graphic  and  tabular  form,  editing,  and  interpolating, 
data  curves.  Digitization  and  Graphics  require  a Tektronix  4394 
or  4395  graphics  tablet  and  4014  Terminal  respectively.  The  other 
tasks  can  be  done  at  any  terminal. 

A complete  description  and  instruction  on  DIGIPLT  is 
given  in  the  DIGIPLT  USERS  GUIDE  in  the  Appendix  G of  this  report. 

In  work  performed  under  this  contract,  two  relatively 
small,  specific-purpose  data  processing  programs  were  written. 

These  programs  are  listed  here  along  with  a brief  description  of 
each. 

DATPLOT  - A program  to  produce  log- log  plots,  with 
axis  systems,  of  crack  growth  rate  versus 
stress  intensity. 

CONVRT  - A program  to  convert  data  output  from  the 
electron  microprobe  to  usable  format;  edit 
erroneous  points;  and  compute  the  mean  and 
standard  deviation  for  each  data  set  pro- 
cessed. 

These  programs  are  listed  in  Appendix  H. 
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5.1.2  LSI-11  Peripheral  Driver  Programs 


The  DEC  RT-11  operating  system  can  be  used  to  generate 
stand  alone  FORTRAN  program  which  contain  a driver  for  a terminal. 

To  drive  other  peripherals  (such  as  A/D  converter,  D/A  converter, 
real  time  clock,  and  modem  interface)  it  is  necessary  to  write 
Assembly  Language  programs  that  can  be  called  from  FORTRAN  pro- 
grams. There  is  a linking  scheme  used  by  RT-11  FORTRAN  for  Assembly 
Language  siabroutines . This  scheme  was  used  in  writing  these  pro- 
grams. The  input  and  output  parameters  of  the  subroutines  and  the 
FOPVK/\N  calls  are  shown  in  Table  V.  Listings  of  these  PDP-11 
Assembly  Language  programs  can  be  found  in  Appendix  I. 

5.1.3  PLli  Support  Programs 

The  PLM  compiler  (Intel  8080  high-level  language) 
originally  used  by  AFML/DOC  was  the  compiler  installed  on  the 
General  Electric  Time-Share  System.  This  compiler  had  two  steps 
of  operation.  The  first  step  had  input  of  an  ASCII  source  PLM 
program  file  and  output  of  an  ASCII  file  that  contained  Assembler 
Language  and  compiler  error  information.  The  second  step  had  as 
input  the  output  file  from  the  first  step  and  output  of  an  ASCII 
file  that  contained  memory  and  symbol  information  as  well  as  the 
memory  laoding  information.  It  was  necessary  to  punch  an  ASCII 
paper  tape  of  the  memory  loading  information  and  convert  that  in- 
formation to  a binary  paper  tape  suitable  for  loading  by  the  paper 

loader  suppl ied  by  Control  Logic.  This  conversion  had  been 
done  previously  on  a Control  Logic  System  with  low-speed  tape 
reaco-t  and  was  very  time  consuming  (45  minutes  for  a 4K  program)  . 

In  addition,  there  was  noise  on  the  telephone  lines  which  made 
pon.-fi''c  a correct  ASCII  tape  very  unlikely.  Many  times  the  tape 
'I’f;  lav'-  to  be  punched  four  oi  five  times  to  get  an  error  free 
♦ n r.  •pp'T’P^.p  pr'iigram  was  written  for  the  GE  system  that  used 
the  ASCII  memory  loading  file  as  input  and  created  a binary  output 
tile  compatible  o'ith  the  paper  tape  loader.  A listing  of  this 
{lOQiitm  is  in  Appenlly  J,  75^0  binary  output  file  was  then  dumped 
tc  the  < . T.voi  ' . . . ^ . r .pi-).  This  process  greatly  improved 


TABLE  V 


CALLING  ASSEMBLY  LANGUAGE  SUBROUTINES 
FROM  RT-11  FORTRAN 


Analog  to  Digital  Converter; 

CALL  ADC  (ICH,  IN) 

where  ICH  = Channel  number  to  be  converted  (0  to  31) 

IN  = Integer  equivalent  of  analog  input  on 
A/D  channel  number  ICH 


Digital  to  Analog  Converter; 

CALL  DAC  (ICH,  lOUT) 

where  ICH  = Digital  to  Analog  converter  number  (1  or  2) 

lOUT  = Integer  equivalent  of  desired  voltage 
on  D/A  channel  number  ICH 


Clock; 


CALL  CLOCK  (IFLG,  NHR, 

NMIN,  NSEC,  NMSEC) 

where  IFLG 

= 

0 for  read  clock, 

1 for  set  clock 

NHR 

= 

Integer  number  of 

hours 

NMIN 

= 

Integer  number  of 

minutes  r 

NSEC 

= 

Integer  number  of 

seconds 

NMSEC 

= 

Integer  number  of 

milliseconds 
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turnaround  time  for  minor  program  changes,  especially  in  the  debug 
procedure.  An  option  was  added  to  print  an  octal  symbol  table  rather 
than  the  standard  hexadecimal  symbol  table. 

Also  a FORTRAN  program  was  written  to  scan  through  the 
output  file  cf  the  first  compiler  step  and  print  errors.  A listing 
of  this  progran’  is  an  Ap'pendix  J.  Previously  it  has  been  necessary 
to  search  this  file  for  errors  line  by  line,  which  was  not  only 
slow  and  difficult,  but  inaccurate.  This  program  would  print  each 
line  with  an  error,  the  error  diagnostic,  and  the  total  number  of 
errors.  Botli  of  these  utility  programs  will  be  adopted  to  either 
the  SEL  86  or  the  CDC  6600  for  future  PLM  use. 

5.2  AFML  USER  EDUCATION 

5.2.1  CDC  6600  Users 

To  determine  methods  AFML  could  employ  to  realize  more 
efficient  utilization  of  allocated  computer  resources,  an  analysis 
of  AFML's  use  of  the  CDC  6600  was  made.  The  results  of  this  analy- 
sis indicated  the  following  would  yield  a significant  increase  in 
efficiency  and  throughout  for  prime  shift  processing: 

1.  Limiting  the  Computer  Resource  Units  available 
per  AFML  prime  shift  job  to  100  CRUS.  This  would  prevent  a small 
number  of  large  jobs  expending  a large  percentage  of  AFML's  allo- 
cated resources,  which  then  results  in  low  throughout  long  turn- 
around AFML  jobs  for  the  remainder  of  the  shift. 

2.  AFML  computer  users  requesting  only  the  resources 
required  for  their  computer  jobs  versus  gross  overestimates,  re- 
su! tirur  in  unnecessary  restrictions  on  computer  resources  available 
to  AFMl.  whi  le  such  jobs  are  in  the  input  queue  or  in  execution. 

}.  More  efficient  use  of  permanent  file  space  by  us- 
i'v;;  a'J  five  available  cycles  per  permanent  file;  and  by  maintain- 
io  ; on  ; etauir-ent  rile  only  information  currently  being  processed. 


To  effect  the  above,  it  was  seen  that  a user  education 
effort  •-.*03  i * crfor  t.  ■ <z\^  d ■'monst -ate  to  AFML  users  the 


r 
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ways  such  practices  caould  be  carried  out  and  how  their  use  would 
benefit  the  individual  user.  A class  was  conducted  during  which 
methods  of  efficient  use  of  the  computer  were  demonstrated  and 
discussed. 

A second  analysis  of  AFML  computer  use  was  made  to 
evaluate  the  effectiveness  of  the  effort.  The  results  showed  an 
increase  of  50%  in  prime  shift  throughout,  and  a significant  in- 
crease in  available  permanent  file  space. 

5.2.2  LSl-11  Users 

A two  day  class  attended  by  14  people  was  held  on  the 
applications  of  the  LSI-11  microcomputer  for  potential  LSI-11  users. 
A wide  variety  of  subjects  were  discussed  ranging  from  very  general 
to  very  specific.  A large  emphasis  was  placed  on  I/O  as  this 
applies  to  all  users. 


APPENDIX  A 

LIQUID  ADSORPTION  EXPERIMENT  PLM  PROGRAM 
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APPENDIX  D 

USERS  MANUAL  FOR  THE  DATA  ACQUISITION  IN  THE  LHMEL  FACILITY 


SECTION  1 
PREVIEW 


The  Laser  Hardened  Materials  Evaluation  Laboratory  (LHMEL) 
facility  includes  a dedicated  Hewlett-Packard  9640A  Multiprogramming 
System  and  a communication  link  to  the  SEL86  computer  system.  The 
OZ  software  system  is  a collection  of  programs  residing  both  on  the 
Hewlett-Packard  system  and  the  SEL86  system  that  acquires  experi- 
mental data  and  performs  real  time  analysis  of  data,  permitting  the 
LHMEL  operator  to  rapidly  alter  experimental  parameter  values  for 
best  effect. 

The  experimental  data  that  are  acquired  are: 

1)  Various  experimental  pre-  and  post-run  operating  parameter  values 
(entered  by  the  operator) . 

2)  The  run  startup  time. 

3)  The  laser  power  (20  Hz  sample  rate) . 

4)  The  air  temperature  (20  Hz  sample  rate) , 

5)  The  static  air  temperature  (20  Hz  sample  rate) . 

6)  The  total  air  pressure  (20  Hz  sample  rate) . 

7)  The  time  the  laser  shutter  opens  and  closes. 

8)  The  time  the  burn-through  wire  in  front  of  the  target  is  cut  by 
the  laser  beam. 

9)  The  time  the  burn-through  wire  behind  the  target  is  cut  by  the 
laser  beam. 

The  principle  component  of  the  Hewlett-Packard  system  is  a 21MX 
computer,  operating  under  the  Hewlett-Packard  supplied  RTE-C  operating 
system.  To  begin  an  experimental  run  the  operator  uses  RTE-C  to 
activate  the  program  OZ  and  enters  the  pre-run  operating  parameter 
values  that  have  changed  since  the  last  experimental  run.  The  opera- 
tor then  has  the  program  OZGO  activated  to  control  the  acquisition 
of  the  experimental  data.  Initialization  of  the  laser  control  se- 
quencer by  the  operator  starts  the  program  data  acquisition  sequence. 


When  the  run  terminates,  the  operator  can  use  RTE-C  to  activate 
the  program  OZPST  to  enter  the  post-test  operating  parameter  values. 

If  desired,  the  operator  can  have  the  program  DADMP  activated  to  in- 
sprr't  the  raw  data  acquired  during  the  run.  If  the  data  are  accept- 
able, the  SEL86  program  can  be  activated.  Witliin  this  program  the 
operator  sends  the  raw  data  over  the  communication  link  to  the  SEL86 
computer  and  initiates  the  execution  of  the  analysis  program  SELOZ 
(on  the  SEL86) . • SELOZ  computes  the  following: 

1)  The  run  time  interval  (time  of  shutter  closure  minus  run  startup 
time) . 

2)  The  shutter  time  interval  (time  of  shutter  closure  minus  time 
shutter  opened) . 

3)  The  burn-through  time  interval  (back  burn-through  wire  time 
minus  front  burn-through  wire  time) . 

4)  The  laser  power  curve  over  the  shutter  time  interval. 

5)  The  laser  power  threshold  time  interval  (the  time  when  the  laser 
power  becomes  less  than,  or  equal  to,  the  laser  power  threshold 
going  down  minus  the  time  when  the  laser  power  becomes  greater 
than,  or  equal  to,  the  laser  power  threshold  going  up) . 

6)  The  maximum,  minimum,  and  average  values  of  laser  power  over  the 
shutter  time  interval. 

7)  The  maximum  and  average  values  of  laser  power  over  the  laser 
power  threshold  time  interval. 

8)  The  total  energy  over  the  laser  power  threshold  time  interval. 

9)  The  diameter  and  area  of  the  beam  on-target. 

10)  The  laser  power  on-target  curve  over  the  shutter  time  interval. 

11)  The  average  power  on- target  over  the  shutter  time  interval. 

12)  The  average  power  density  on-target  over  the  shutter  time  interval. 
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i.  ni .ur n jver  the  shutter  time  interval. 

RMS  value  of  static  air  pressure  over  the  shutter  time  interval 
tota’  a’’-  ,-fT-vo  over  the  shutter  time  interval. 

> V j . uw  MJ  ..  -l;  c.  er  f'e  shutter  time  interval. 


The  data  and  analytical  results  are  conveniently  formatted 
and  printed  on  a high-speed  printer  at  the  SEL86  site.  They  are  also 
archived  on  magnetic  tape.  Programs  are  provided  to  produce  CALCOMP 
quality  plots  of  the  power  on-target,  total  temperature,  static 
pressure,  and  total  pressure  curves  either  in  real-time  or  from  the 
archived  data.  A summary  of  the  analysis  is  printed  for  the  operator 
on  the  LHMEL  teletype.  The  operator  can  use  this  summary  to  select 
the  experimental  conditions  for  the  next  run. 
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SECTION  2 
PROGRAM  SUMMARIES 


2 . 1 INTRODUCTION 

There  are  five  sets  of  routines  within  the  OZ  system: 

1)  RTE-C : the  operating  system  on  the  21MX  computer. 

2)  The  OZ  programs:  tne  set  of  21MX  resident  routines  that  acquire 
the  pre-run  and  post-run  operating  parameter  values  and  experi- 
mental data. 

3)  SEL86 : the  21MX  resident  routines  that  communicate  with  the 
SEL86  computer,  sending  data  to  it,  starting  programs  on  it,  and 
receiving  ASCII  data  from  it  for  output  on  the  LHMEL  teletype. 

4)  SELOZ:  the  set  of  data  analysis  programs  resident  on  the  SEL86 . 

5)  TPOUT:  the  program  resident  on  the  SEL86  that  extracts  data 
archived  on  magnetic  tape  and  offers  a plotting  option. 

Each  of  these  sets  of  routines  will  be  briefly  described  in 
the  following  sections.  A more  detailed  description  is  found  in 
Reference  15., 

2 . 2 RTE-C 

RTE-C  is -Hewlett-Packard's  acronym  for  Real-Time  Core-Based 
Software  System.  The  procedure  for  initializing  this  system  is 
found  in  Section  3.2.  Once  RTE-C  is  running,  and  any  time  it  has 
control  of  the  teletype,  entering  any  character  through  the  tele- 
type v-all  cause  RTE-C  to  respond  by  typing  an  asterisk  (*)  . At  this 
the  operator  is  free  to  enter  any  command  described  in  the 
PTK-C  operating  manual.* 

The  operator  is  required  to  enter  the  time  of  day  every  time 
iy  ’niri^lized  The  format  for  this  command  is: 

TM , day , h , mi n , sec 

-»  i rhree  -di  ait  '^av-of- the-year  (see  Table  2-5  of  the 
i.'o.-atrng  manual),  and  h, min, sec  is  the  current  time  on  a 24- 


I .♦  \ • 


■ hvstem  Programming  and 

T ■ nDOl . 


Among  the  many  functions  of  RTE-C  is  the  control  of  the  opera- 
tor's access  to  the  programs  stored  in  the  21MX  core.  In  order  to 
execute  a program,  the  operator  types: 

ON , pgnam 

t 

where  pgnam  is  the  name  of  the  program  to  be  executed. 

Error  messages  from  RTE-C  are  discussed  on  page  3-12  of  the 
operating  manual. 

It  should  be  noted  that  if  a wrong  character  is  typed  "control- 
A"  can  be  used  to  delete  that  character.  RTE-C  types  a left-facing 
arrow  to  indicate  that  the  character  has  been  deleted.  To  delete 
an  entire  line,  type  "rubout".  RTE-C  types  a back  slash  to  indicate 
that  the  line  has  been  deleted. 

2.3  THE  OZ  ROUTINES 

2.3.1  Introduction 

The  following  OZ  routines  can  be  called  using  the  RTE-C 
ON, pgnam  command: 

OZ 

OZGO 

DADMP 

OZPST. 

In  addition,  two  other  routines  are  called  internally  by  RTE-C: 

OZAD 
OZUNI . 

2.3.2  OZ 

The  OZ  routine  permits  entry  of  pre-run  operating 
parameter  values.  When  the  program  is  started  by  RTE-C,  this  head- 
ing is  typed  on  the  LHMEL  teletype: 

OZ,  VERSION  xx/xx/xx. 

After  this,  the  system  waits  for  the  operator  to  enter  a two-letter 
code  corresponding  to  a pre-run  operating  parameter.  (These  codes 
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are  listed  in  Table  I.)  The  code  should  be  followed  by  a carriage 
return. 

I 

OZ  searches  its  internal  table  of  codes.  If  it  finds 
a match  it  types: 

ENTER  XX 

where  XX  is  an  echo  of  the  operator- supplied  code.  A carriage  re- 
turn is  executed  and  OZ  waits  for  the  operator  to  enter  the  appro- 
priate parameter  value  followed  by  a carriage  return. 

If  no  match  is  found,  OZ  types: 

NO  MATCH 

followed  by  a carriage  return  and  waits  for  the  code  to  be  reentered. 

Any  number  of  parameters  and  their  values  can  be  enter- 
ed during  an  OZ  execution,  and  their  order  is  arbitrary.  This  per- 
mits entry  of  only  those  parameters  that  have  changed  value  since 
the  previous  run. 

If  a parameter  value  is  entered  incorrectly,  the  param- 
eter code  and  value  can  be  reentered. 

All  numeric  parameter  values  should  be  entered  in  F7.3 
format  (i.e.,  in  the  form  999.999).  For  alphanumeric  parameter 
values,  the  length  of  the  field  varies  with  parameters  as  defined 
in  Table  I. 

To  notify  OZ  that  no  more  parameters  are  to  be  entered, 
type  XX  for  the  parameter  code.  OZ  types: 

END  OZ 

r‘ tar  nr  ‘ control  to  RTC-C. 

I'lC  CZ,0  initiates  the  data  acquisition  process. 

■ •iL'?ri  .q  *-be  routine,  all  the  flags  and  pointers  are  initialized 
and  the  rear  bujn-thiough  wire  clock  is  zeroed.*  The  following  mes- 
J.S  printed  o.i  the  Loir  type. ; 

^poir'rVrmi' .'.yt  1-.  1 ot  a r.  ■ t 
tha!;  t’  . ' ■ 


r'iOT  oa\n-througn  wire  time  indicates 


TABLE  I 

PRE-RUN  OPERATING  PARAMETERS 


PARAMETER 

CODE 

FORMAT 

UNITS 

sample  id 

SI 

A6 

test  number 

TN 

A4 

. 

distance  between  mirror  and 

sample  DI 

F7.3 

cm 

HE  dome  pressure 

HD 

F7.3 

psig 

N2  dome  pressure 

ND 

F7.3 

psig 

C02  dome  pressure 

CD 

F7.3 

psig 

laser  pressure 

LA 

F7.3 

. mm  Hg 

exhaust  pressure 

EP 

F7.3 

mm  Hg 

high  voltage 

KV 

F7.3 

kilovolts 

run  time 

RT 

F7.3 

sec 

HE  line  pressure 

HL 

F7.3 

psig 

N2  line  pressure 

NL 

F7. 3 

psig 

C02  line  pressure 

CL 

F7.3 

psig 

WT  back  pressure 

WB 

F7.3 

psig 

WT  line  pressure 

WL 

F7.3 

psig 

WT  plenum  temperature 

WP 

F7.3 

deg 

WT  nozzle  design,  x dim 

WX 

F7.3 

cm 

WT  nozzle  design,  y dim 

WY 

F7.3 

cm 

initial  load 

IL 

F7.3 

kilopounds 

I 

I 

I 
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READY. .AUTO? 


The  operator  should  respond  by  entering  a "Y"  or  an  "N"  followed 
by  a carriage  return. 

If  a "Y"  is  entered,  the  system  waits  up  to  20-seconds 
for  a signal  from  the  laser  control  system  (on  channel  0 of  the 
universal  interface  card  data  buffer) . 

(a)  If  the  pulse  does  not  occur  within  20-seconds, 
this  message  is  typed  on  the  teletype; 

TIME  UP 

and  control  is  returned  to  RTE-C. 

(b)  If  the  pulse  does  occur  within  the  time  limit, 
the  OZUNI  routine  sets  the  experiment  initialize  flag  IF4  and  re- 
cords the  time  of  the  pulse  as  the  "start-up"  time.  Control  is  re- 
turned to  RTE-C. 

If  an  "N"  is  entered,  the  experiment  initialize  flag 
IF4  is  set,  the  current  time  is  recorded  as  the  "start-up"  time, 
and  OZAD  and  OZUNI  are  started.  Control  is  returned  to  RTE-C. 

2.3.4  DADMP 

The  DADMP  routine  generates  a teletype  listing  of  the 
values  of  pre-run  and  post-run  parameters,  the  flags,  the  pointers, 
the  data  arrays  and  the  time  buffers. 


The  numbers  under  the  "FLAGS  AND  POINTERS"  banner  are 

lJiu  /c  j ufc.j  ol  the  flags  and  pointers  in  this  sequence; 

y . vpnv.i.ir.crt  initialize  flag;  incremented  each  time  the  experi- 
mert  is  Lritial.ized  during  a run. 

i’ ‘ bu!  '.  tl  t 'ugh  wi’-e  flag,  incremented  each  time  a front 
i tnrovjh  'jfir.  rijral  is  detected  during  a run. 

- •<  i'  bu'-n- Through  wire  flag;  incremented  each  time  a rear 

-j'.iij ' i-AiiOugii  witv  bJ.gjial  is  detected  during  a run. 

'.butter  ocen  flaa;  incremented  each  time  a shutter  open  signal 
is  detected  during  a run. 


IFER 


- previous  universal  interface  data  register  contents  (decimal 
equivalent) . 

IRBUF-  current  universal  interface  data  register  contents  (octal 
value) . 

IFIN  - set  when  flags  initialized  and  reset  when  run  completed. 

IPNTR-  the  number  of  filled  entries  in  the  data  arrays  (decimal 
value) . 

The  sequence  of  numbers  labeled  "START  UP="  down  through 
"SHUTTER  CLOSE="  are  associated  with  the  clock  buffer  entries.  The 
order  of  these  entries  on  a single  line  is:  tens  of  msec,  seconds, 
minutes,  hours,  days  (all  in  decimal) . 

The  values  in  the  data  arrays  are  in  units  of  volts. 

2.3.5  OZPST 

The  OZPST  routine  is  used  to  enter  the  post-run  opera- 
ting parameter  values.  When  the  program  is  started  by  RTE-C,  the 
following  heading  is  printed  on  the  LHMEL  teletype: 

OZPST, VERSION  XX/XX/XX. 

At  this  point  the  system  idles,  waiting  for  the  operator  to  enter 
one  of  the  two-letter  codes  corresponding  to  a post-run  operating 
parameter  listed  in  Table  II.  The  code  should  be  followed  by  a 
carriage  return.  ' 

OZPST  searches  its  internal  table  of  codes.  If  a match 
is  found,  it  types 

ENTER  XX  ' 

r 

where  XX  is  the  echo  of  the  operator-supplied  code.  A carriage  re- 
turn is  executed  and  OZPST  waits  for  the  operator  to  enter  the 
appropriate  parameter  value  followed  by  a carriage  return. 

If  no  match  is  found,  OZPST  types 

NO  MATCH 

followed  by  a carriage  return  and  waits  for  the  code  to  be  reentered. 

r 

All  the  parauneter  values  are  numeric  and  should  be 
entered  in  an  F7.3  format  (i.e.,  in  the  form  999-999). 
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TABLE  II 

POST-RUN  OPERATING  PARAMETERS 


PARAMETER 

CODE 

FORMAT 

UNITS 

HE  upstream  pressure 

HU 

F7.3 

psig 

N2  upstream  pressure 

NU 

F7.3 

psig 

C02  upstream  pressure 

CU 

F7.3 

psig 

laser  pressure 

LP 

F7.3 

mm  Hg 

exhaust  pressure 

EX 

F7.3 

mm  Hg 

current 

CR 

F7.3 

mi  11 lamps 

WT  exhaust  pressure 

WE 

F7.3 

psig 

i: 

i 


1 


r 


To  notify  OZPST  that  no  more  parameters  are  to  be 
entered,  the  operator  should  type  "XX"  for  the  parameter  code. 

OZPST  types: 

END  OZPST 

and  returns  control  to  RTE-C. 

2.3.6  OZAD 

The  purpose  of  OZAD  is  to  control  the  A/D  converter. 

When  OZAD  is  started  by  OZUNI,  it  checks  the  IF4  flag 
(experiment  initialize)  to  determine  if  the  run  has  been  initialized 
by  the  OZGO  routine.  If  it  has  not,  this  message  is  typed  on  the 
teletype : 


OZAD  NOT  INIT. 

Control  is  returned  to  RTE-C. 


If  the  run  has  been  initialized,  the  IF3  flag  (shutter 
closed)  is  examined  to  determine  if  the  run  has  been  cpmpleted.  If 
it  has,  the  message 

A/D  DONE  : 


is  typed  on  the  teletype  and  control  is  returned  to  RTE-C.  If  the 
run  is  still  in  progress,  a single  reading  is  taken  on  each  of  the 
first  four  channels  (CH0-CH3)  of  the  A/D  converter.  These  readings 
are  stored  in  the  next  available  location  of  the  four  arrays 
PMETR(CHO),  TTEMP(CHl),  SPRES(CH2),  and  TPRES(CH3).  Then  the  pointer 
used  to  indicate  the  next  available  location,  IPNTR,  is  incremented 
and  that  value  compared  with  the  maximum  dimension  of  the  four  arrays 
(101) . If  the  arrays  are  not  full,  OZAD  reschedules  itself  to  exe- 
cute 50  msec  later  and  control  is  returned  to  RTE-C. 

If  the  arrays  are  full,  which  occurs  five  seconds  after 
the  run  begins,  this  message  is  printed  on  the  teletype: 


ARRAYS  FULL 
A/D  DONE 

Control  is  returned  to  RTE-C. 


1 
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2.3.7  OZUNI 


OZUNI  is  the  program  that  processes  the  signals  input 
to  the  Universal  Interface  Card  (12930A) . It  is  started  by  the 
OZGC  program  anvl  it  reschedules  itself  to  run  every  20  msec  until 
either  the  ran  is  completed,  an  error  has  occurred,  or  20  seconds 
has  passed  wxthout  a "start-up"  signal.  This  program  should  not 
be  called  by  the  operator  using  the  RTE-C  "ON, OZUNI"  command. 

The  Universal  Interface  Card  data  channel  input  register 
is  assumed  by  OZUNI  to  be  wired  as  follows: 

bit  0:  run  start-up  signal  wire 
bit  1:  shutter  signal  wire 
bit  2:  front  burn-through  signal  wire 
bit  3:  back  burn-through  signal  wire. 

All  signals  are  assumed  to  be  single-ended- 

The  logical  meaning  of  each  of  the  signals  is  assumed 
to  be  as  follows; 


Run  start-up  signal:  the  occurrence  of  the  0 to 
5v  transition  indicates  start  of  run.  The  signal 
remains  at  5v  for  the  duration  of  the  run. 

Shutter  signal  wire;  the  occurrence  of  the  first 
0 to  5v  transition  after  the  start  of  a run  in- 
dicates the  opening  of  the  shutter.  Occurrence 
of  the  next  5v  to  Ov  transition  indicates  the 
closing  of  the  shutter.  The  signal  remains  at 
5v  in  between. 

Front  burn-through  signal:  the  occurrence  of 
the  first  0 to  5v  transition  after  the  start  up 
of  a run  indicates  that  the  front  burn-through 
v.a.re  has  burned  through.  The  signal  remains  at 
b'j  it  I the  duration  of  the  run. 

Back  burn-through  signal:  the  occurrence  of  the 
first  0 to  Sv  transition  after  the  start  of  a run 
i.i.  ii  ates  hat  the  back  burn-through  wire  has  burn- 
ed -.lirougl).  The  signal  remains  at  5v  for  the  dura- 
tion of  Hie  run. 

Upon  the  first  entry  into  the  program  during  a run, 
tni;  i'lNIT  flag  is  examined  to  determine  if  the  OZGO  program  ini- 


EARLY  UNIV  INT 


is  typed  on  the  teletype  and  control  is  returned  to  RTE-C.  This 
I prevents  the  operator  from  starting  OZUNI  using  the  RTE-C  "ON,OZUNI" 

S command.  Otherwise,  each  time  OZUNI  runs,  the  first  thing  it  does 

is  load  the  Universal  Interface  Card's  data  register  into  the  A 
li  register,  where  it  is  exclusive-ored  with  the  contents  of  the  data 

register  obtained  the  last  time  OZUNI  ran.  The  resulting  word  has 
a bit  set  in  every  position  where  the  level  of  the  data  register 
has  changed  since  the  last  time  OZUNI  ran.  This  word  is  processed 
as  shown  in  the  accompanying  flowchart  in  Figure  1.  Djuring  the 
course  of  processing,  if  a logical  error  occurs,  (e.g.,,  two-level 
I changes  in  the  front  burn-through  wire  signal  during  t^e  course 

1 of  one  run)  , the  message  "ERROR  IN  UNIV  INT"  is  printe,d  on  the  tele- 

I type  and  data  acquisition  for  the  run  is  aborted.  The  source  of 

1 the  error  can  be  determined  by  running  DADMP  and  examining  the 

values  for  IFER  and  IRBUF.  If  no  logical  error  occurs,  OZUNI 
checks  the  shutter-closed  flag  to  determine  if  the  shutter-closed 
event  has  occurred.  If  it  has,  the  program  terminates.  If  it 
hasn't  the  program  reschedules  itself  to  run  20  msec  later. 

t 

After  OZGO  turns  OZUNI  on,  OZUNI  runs  ev^ry  20  msec, 

: checking  for  the  start-up  signal.  OZGO  times  the  evei^t.  If  the 

j signal  is  not  generated  in  20  seconds,  OZGO  notifies  OZUNI  through 

the  IFIN  flag,  and  OZUNI  does  not  reschedule  itself. 

I 

2.4  SEL86 

• The  SEL86  routine  controls  the  user's  communication  with  the 

SEL  86  computer.  SEL86  is  brought  up  by  typing. 

ON,SEL86. 

I The  response  of  the  terminal  is  ‘ 

-SEL  TERMINAL  SUPPORT  SYSTEM,  TEIiMINAL  TY00  - 
ENTER  FUNCTION  CODE,  ?,  OR  CR  TO  TERMINATE; 

The  informed  user  will  recognize  this  as  the  normal  mode  of  entry 
into  the  SEL  86  Terminal  Support  System  (TSS) . All  the  options 
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PROCESS 

START-UP 

SIGNAL 


PROCESS 

SHUTTER 

SIGNAL 


t 1.  B 


i \ c. 


\ 


of  this  system  are  now  open  to  the  user.  The  normal  protocol  for 
an  experiment  is  to  use  TSS  to  transfer  the  COMMON  area  of  RTE-C 
to  a disk  file  on  the  SEL  86  and  then  to  run  SELOZ. 

To  transfer  the  COMMON  area,  type  "RM"  in  response  to  the  above 
request: 

ENTER  FUNCTION  CODE,?,  OR  CR  TO  TERMINATE:  RM 

The  terminal  will  respond  with  a question  mark  and  the  user  should 
type  Y:  ' 

?Y 

The  terminal  will  follow  this  with  another 

ENTER  FUNCTION  CODE,?,  OR  CR  TO  TERMINATE: 

The  user  should  respond  with: 

ENTER  FUNCTION  CODE,?,  OR  CR  TO  TERMINATE:  PM. 

The  terminal  will  type: 

- BEGIN  PROGRAM  MONITOR  - 
PROGRAM  (,  # SECONDS ) : 

The  proper  user's  response  is: 

PROGRAM  (,  #SECONDS) : M0SELOZ, 200. 

This  will  start  the  execution  of  SELOZ. 

When  SELOZ  has  completed  execution,  the  normal  sequence  of  inter- 
action is 

PROGRAM (, #SECONDS) : control-D 

ENTER  FUNCTION  CODE,?,  OR  CR  TO  TERMINATE:  control-D 
TERMINAL  TY00  LOGGED  OUT 
At  this  point  control  is  returned  to  RTE-C. 

At  the  end  of  a day's  runs,  provide  the  following  command  to 
RTE-C : 

*OF,SEL86,l. 

A more  technical  description  of  SEL86  is  found  in  the  Pro- 
grammer Comments  section  of  Reference  15. 
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2.5  SELOZ 


SELOZ  and  its  supporting  routines  analyze  the  data  obtained 
by  OZ.  generate  summaries  of  the  analysis  produce  plots,  and  store 
the  data  on  digital  magnetic  tape. 

Wiien  SELOZ  is  initiated  through  SEL86  it  performs  the  follow- 
ing procedure: 

1)  Send  the  message  "BEGINNING  SELOZ"  to  the  LHMEL  teletype. 

2)  Input  the  Calibration  Deck  (M0$CALDK) . 

3)  Input  the  experimental  data  (M0$OZFIL) . 

4)  Calculate  the  time  the  shutter  was  open  (TS) : 

TS  = shutter  close  time  (ISHC)  - shutter  open  time  (ISHO) . 

5)  Calculate  the  time  required  to  burn  through  the  sample  (TB) ; 

TB  = back  burn  through  wire  time  (IBTB)  - 
front  burn-through  wire  time  (IBTF) . 

6)  Calculate  the  time  between  the  shutter  open  time  and  the  time 
the  front  burn-through  wire  burned  through  (TD) ; 

TD  = front  burn-through  wire  time  (IBTF)  - 
shutter  open  time  (ISHO) . 

7)  If  TD  is  greater  than  30  seconds  or  negative,  dump  results, 
which  will  consist  of  zero  values,  and  exit. 

8)  Convert  the  power  meter  readings  to  floating  point  volts: 

CH0(I)  = FLOAT(IAND(IMETR(I) , B177760) )*. 0003125  . 

(Note:  This  procedure  is  documented  in  the  HP  91000A  A/D 

Convex  ter  Manual . ) 

9^  Convert  power  meter  volts  to  laser  power  (PL(I)): 

PL(J)  = power  volts  + (power  responsivity /beamsplitter 
reflectivity) 


or 

rr.ii;  - ^ (c(2)/c(l)) 

( i.;e  T;’ble  TIT  for  Calibration  Table  parameter  definitions  as 
j - -ill  r-r  ’ ;3  arrays.) 

Convert  laser  power  to  laser  power  on  target  (PT(I)); 

: C r-.:.f  1 PCti  vi  ty)  * 

Li  / L . j ...f  r-.irrors)  * 

I. 


10> 


or 

PT(I)  = {l.-C(L))*  C(3)  * PL(I) 

11)  Find  the  time  (TLl)  at  which  laser  power  (PL)  first  equals  or 
exceeds  the  laser  power  threshold  (C(4)). 

12)  Find  the  first  time  (TL2)  after  TLl  at  which  the  laser  power 
(PL)  decreases  below  the  laser  power  threshold  (C(4)). 

13)  Calculate  the  difference  between  TLl  and  TL2  (TL) : 


TL  = TL2-TL1  . 


14)  Define  the  TL  time  interval  as  the  time  interval  from  TLl  to 
TL2  inclusive. 

IBINLl  is  the  element  in  any  of  the  four  data  arrays  correspond- 
ing to  TLl  and  IBINL2  is  the  element  corresponding  to  TL2.  Then 
NL  = IBINL2-IBINL1+1  is  the  number  of  elements  in  any  array  over 
the  TL  time  interval.  Calculate  the  average  laser  power  over 
the  TL  time  interval  (PLLA) ; 


, IBINL2 

PLLA  = ^ Z PL(I)  . 

I=IBINL1 


15)  Calculate  the  maximum  laser  power  (PLLMX)  over  the  TL  time 
interval ; 

PLLMX  = max  {PL(I)  I I=IBINL1,  IBINL2}  . 

16)  Define  the  TS  time  interval  as  the  time  interval  from  shutter 
open  time  (ISHO)  to  shutter  close  time  (ISHC)  inclusive.  De- 
fine IBINSl  as  the  element  in  any  of  the  four  data  arrays 
corresponding  to  ISHO  and  IBINS2  as  the  element  corresponding 
to  ISHC.  NS=IBINS2-1BINS1+1  is  the  number  of  elements  in  any 
data  array  over  the  TS  time  interval.  Calculate  the  average 
laser  power  over  the  TS  time  interval  (PLSA) : 


PLSA 


1 

NS 


IBINS2 

Z PL(I)  - 

I=IBINS1 
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17)  Calculai.e  the  maximum  laser  power  over  the  TS  time  interval 
(PLSMX) : 

max  {PLfl ) ! 1=IBINS1,  IBINS2}  . 

18)  Calculate  the  minimum  laser  power  over  the  TS  time  interval 
(;',L.SMN)  : 

PLSMN  = min{PL(I) I I=IBINS1,  IBINS2}  . 

19)  Calculate  the  integral  of  laser  power  over  the  TL  time  inter- 
val (WJ)  : 

IBINL2 

WJ  = .050*  I PL(I)  . 

I=IBINL1+1 

The  value  .050  is  in  units  of  seconds  and  is  the  bin  width. 

20)  Calculate  the  average  laser  power  on  target  over  the  TS  time 
interval  (PTSA) ; 


PTSA  = 


NS 


IBINS2 

Z PT(I) 

I=IBINS1 


21)  Define  the  TB  time  interval  as  the  interval  from  the  time  the 
front  burn-through  wire  burned  through  to  the  time  the  back 
burn  through  wire  burned  through.  Define  IBINBl  as  the  element 
in  any  of  the  four  data  arrays  corresponding  to  IBTF  and  IBINB2 
as  the  element  corresponding  to  IBTB.  NB  - IBINB2-IBINB1  is 
tae  number  of  elements  in  any  data  array  over  the  TB  time  in- 
ti.aval.  Calculate  the  average  laser  power  on  target  during 
t’  ^ TP  time  interval  (PTBA)  ; 


’TB  A “ 


1 

Nfl 


IBINB2 

Z PT{I) 

I=--IBINB1 


C > ipr  ,-c  ‘he  mirror-to-t'3mple  length  supplied  in  the  experi- 
.’i.untal  data  (FID(l))  with  the  three  sets  of  distance  bounds 
in  tnc  I •:  io'  Tab'e.  Select  the  correct  slope  (C7)  and 

t L ; ■ . ' . . . - ..e  M,';  area  of  the  laser  beam 


t 


at  the  target: 

DIAM  = C6+C7*FID(1) 

AREA  = (j)*DIAM^  . 

23)  Calculate  the  average  laser  power  on  target  per  unit  area  of 
the  target  over  the  TS  time  interval  (PDSA) : 

PDSA  = PTSA/AREA  . 

24)  Calculate  the  average  laser  power  on  target  per  unit  area  of 
the  target  over  the  TB  time  interval: 

PDBA  = PTBA/AREA  . 

25)  Calculate  the  integral  of  laser  power  on  target  over  the  TS 
time  interval  (EDSA) : 

IBINS2 

EDSA  = .050*  I PT(I) . 

I=IBINS1+1 

f 

26)  Calculate  the  integral  of  laser  power  on  target  over  the  TB 
time  interval  (EDBA) : 

IBINB2 

EDBA  = .050*  I PT(I)  . 

I+IBINBl+1 

27)  Convert  the  air  temperature  readings  to  floating  point  volts 

CHl(I)  = FLOAT(IAND(ITEMP(I) , B177760) )* .0003125  . 

28)  Convert  air  temperature  volts  to  degrees 

air  temperature,  degrees  = (air  temperature,  volts  ) + 
air  temperature  channel  responsivity,  degrees/volt) 

or 

CT(I)  = CH1(I)*C(12)  . 

29)  Convert  the  static  pressure  readings  to  floating  point  volts 


CH2(I)  = FLOAT ( I AND (I PRES ( I) , B177760) )*  .0003125 


30)  Convert  static  pressure  volts  to  psi: 


static  pressure,  psi  = (static  pressure,  volts)* 
(static  pressure  channel  responsivity,  psi/volt) 


or 


CT2(I)  = CH2(I) *C(13)  . 

31)  Con\ert  the  total  pressure  readings  to  floating  point  volts: 

CHJ(I)  = FLOAT(IAND(JPRES(I) ,B177760) ) *.0003125  . 

32)  Convert  total  pressure  volts  to  psi: 

total  pressure,  psi  = (total  pressure,  volts)* 

(total  pressure  channel  responsivity,  psi/volt) 

or 

CT3(I)  - CH3(I)*C(14)  . 

33)  Compute  the  minimum  air  temperature  over  the  TS  time  interval 
(CTMN) : 

CTMN  = Mih  {CT(I)  !I=IBINS1,  IBINS2}  . 

34)  Compute  the  RMS  value  of  the  static  pressure  over  the  TS  time 
interval  (PSR) : 


PSR 


- 

y 1= 


1BIH52 

I CT2(I) 
IBINSl 


■)  C’.lculate  the  RMS  value  of  total  pressure  over  the  TS  time 
interval  (PTR) : 


PTR 


- 

V NS 


IBINS2 

I CT3(I) 

I=IBINS1 


< . l..  i .1  n he  rac/i  numt  (CMACH)  : 


CMACP  - - 1.1 


37)  Print  the  Calibration  Table. 

38)  Print  the  evaluation  sheets  and  dumps  of  various  arrays.. 

39)  Send  the  post- test  summary  to  the  LHMEL  teletype. 

40)  Send  "SAVE  DATA?"  message  to  LHMEL  teletype.  Wait  for  an  answer. 
If  the  answer  is  not  "N",  sent  the  message  "NEW  TAPE?"  to  the 
LHMEL  teletype.  If  the  response  to  "NEW  TAPE?"  is  "N",  send 

the  response  to  "FIRST  RUN?"  to  the  LHMEL  teletype.  If  the 
response  to  "FIRST  RUN?"  is  "N",  then  this  is  not  the  first  run 
of  the  day  and  the  run  data  are  added  to  the  tape.  If  the 
response  to  "FIRST  RUN?"  is  anything  else,  then  a search  of  the 
magnetic  tape  is  made  for  the  last  run  on  the, tape  before  the 
run  data  are  added  to  the  tape.  If  the  response  to  "NEW  TAPE?" 
is  anything  other  than  "N"  then  the  tape  is  assumed  to  be  empty 
and  the  run  data  is  written  without  chec)cing  whether  or  not  the 
tape  is  at  the  end  of  useful  information. 

41)  Send  "PLOTS?"  message  to  LHMEL  teletype.  Wait  for  an  answer. 

If  the  answer  is  not  "N",  generate  the  plots. 


i 


TABLE  III 

CT^LIBRATION  TABLE  PARAMETERS  IN  ARRAYS 


ARRAY  LABEL  CONTENTS 

DATE  calibration  table  date 

C(l)  beamsplitter  reflectivity 

C(2)  power  responsivity  (watts/v) 

C(3)  combined  reflectivity  of  3 mirrors 

C(4)  laser  power  threshold  (watts 

C(5)  mirror  focal  length  (m) 


Did) 

D2(l) 

C(6) 

C(7) 

Dl(2) 

D2(2) 

C(8) 

C(9) 

Dl{3) 

D2(3) 

C(10) 

C d 1 ) 

Cd2) 

Cvi3) 

■-xx4) 


lower  mirror  distance  bound  'j  diameter  parameters  used 

, , (cm)  > for  sample  palced  much 

upper  mirror  distance  bound  j ^han  focal  lengths 

slope 

intercept 

lower  mirror  distance  bound  'j  diameter  parameters  used 

(cm)  > for  sample  placed  near 
upper  mirror  distance  bound  J focal  lengths 

slope 

intercept 

lower  mirror  distance  bound  , diameter  parameters  used 

(cm)  J for  sample  placed  much 

upper  mirror  distance  bound  J greater  than  focal  lengths 

slope 

intercept 

air  temperature  channel  responsivity  (degree/volt) 
static  pressure  channel  responsivity  (psi/volt) 
total  pressure  channel  responsivity  (psi/volt) 


SECTION  3 

STANDARD  OPERATING  PROCEDURES 

3.1  BASIC  PROCEDURE 

1.  Before  each  run  in  which  the  tape  drive  or  plotter  is  to 
be  used,  or  before  each  series  of  runs  if  the  runs  are  going  to 
follow  one  after  another,  call  the  SEL  operator  (Cindy)  at  53778 
and  inform  her  of  the  equipment  requirements.  Also  inform  her 
whether  or  not  the  tape  is  to  be  rewound  at  the  end  of  the  run. 

2.  Start  RTE-C  using  Section  3.2. 

3.  Load  the  OZUNI  program  using  the  RTE-C  "LO"  command. 

4.  If  necessary,  modify  the  Calibration  Table  using  the 

procedure  in  Section  3.3.  ; 

5.  Bring  up  OZ  to  set  the  pre-run  parameters  values. 

6.  Bring  up  OZGO  for  the  experimental  run. 

7.  Bring  up  OZPST  to  set  the  post-run  parameter  values. 

8.  If  desired,  bring  up  DADMP  to  examine  the  data. 

9.  Bring  up  SEL86  and  transfer  the  data  to  the  OZFIL  file 
in  the  SEL. 

10.  Run  SELOZ. 

11.  Bring  up  OZ  to  change  any  pre-run  parameters  that  are  to 
be  changed  for  the  next  run . 

12.  Go  to  step  6 for  another  run. 

3 . 2 STARTING  RTE-C 

1.  Press  Master  Power  button  on  upper  right  side  of  9640A 

rack . 

2.  Insert  Key  into  front  panel  of  21MX.  Turn  key  counter- 
clockwise to  R and  then  back  to  operate. 

3.  Turn  on  paper  tape  reader. 

4.  Turn  on  teletype  (to  line  position). 

5.  Load  RTE-C  paper  tape  into  paper  tape  reader. 

6.  Enter  OOlSOOg  into  switch  (5)  register. 

7.  Press  IBL. 
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8.  Press  RESET. 

9.  Press  RUN. 

10.  When  the  paper  tape  has  been  read,  the  light  pattern  on 
rti.at  parcel  ohuaid  Lu  100277  and  this  should  be  the  contents 


't  of  the  S register. 

1 1.1  Enter  000002g  inco  P register. 

12.  Press  Preset  and  RUN. 


3.3  MODIFYING  THE  CALIBRATION  TABLE 

The  Calibration  Deck  used  by  the  SELOZ  program  can  be  changed 
from  the  LHMEL  teletype  using  the  SEL86  TEST  EDITOR.  Any  time  after 
RTE-C  has  been  started , an 

ON,SEL86 

command  will  initialize  communication  with  the  SEL.  The  following 
sequence  will  provide  access  to  the  M0$CALDK  Edit  file,  which  can 
then  be  changed  using  the  EDITOR  commands: 

-SEL  TERI«!1NAL  SUPPORT  SYSTEM,  TERMINAL  TY00 
ENTER  FUNCTION  CODE,?,  OR  CR  TO  TERMINATE:  TO 

. -BEGIN  TEST  EDITOR. 

ENTER  2-CHARACTER  USER  CODE  OR  CR  TO  TERMINATE:  M0 

COMMAND?  USE  CALDK 

COMMAND? 

V.hjo  the  change-s  have  been  made,  the  following  response  saves  the 
file: 

COMMAND?  SAVE  CALDK  SCR 

H .estcn;--'  r-f  E”IT  terminates  the  EDITOR. 

Ic'-i’rj  .:firirra'.d  termi.nates  SEL86: 

OF,SCr..86,l  . 


••..I.l*."  I I •'  I'l  I ‘ 

"O  js'.l**! 

•lilNl  ?«><•»  IHi.li » , 1 s ‘ j • > j rtf  I f (1 ) » IHK ( J ) 


i 

i 

1, 

I 

p 

1 

i 

i 


103 


su'^PouTi’^F  of’Tst  FTN  nf/oi/T7  paof 


tp«;t  ruFM^MT 


J 


«■!)  ifj;  ijpnsrf  ( r . j.  It . f y.  T r I irov,’’) 

01  '-'ISTON  Tin  .'’I  l) 

I'iTTr.ro  n,T 


Slin^o  ITIM'  A . '•>  . TI  n.  V ) 


u 

o 

o 


SUl!?0Uri'4''  7i*/’4  0PT  = 1 FTN  4. ^♦414  . O'?. 74. 51  PSG*’ 


J- )1 
TT=TT-1 

catv  510'??  (n.n.c.GTiCK.T . j.n.  t’')h.'’pt'jT7) 


7J*/7U  T^T  = 1 '4.'«‘t*l4  0f»/'^l/77  0'?*’4.51  o/VGF 


SHIFT  P^COCOS  IN  C0®E!TMONT  HFTT  '^rrQon 


SU^?OUTIN':  TESTtn,n,C,GrAnK,TtJP.L  ,NP0WS,NW^.IY.TY1  , 
dtm-nston  0(1)  .'itH.Cd)  .oTArYiiiinn.ji 

TNTPGF^  0,O,r.STflr.K  , rnr.PnitjT-* 


;ST  UO'(FP  RIGHT  ELEH'’nT 


poor.oftM  7U/7U  rT’4  r^/'^i/77  oflGP 


I ^ I I I 
sD  ^ 


■ y*-ft^are  t t t r t : 
: u zcz^-oo'o  Cl  O'  ct  \ 
ckrc^-ca  »-ooocc< 

i c:(r-t^Q.i/^iiuu.u.ui 


U*?*"  TNTl?.^  TO  SFA.^CH  AMO  ( CO^vr  f^«;X  ^N) 


►-  c 

V u.- 

c L r- 


u t I u 
T O 0 
ll  « 

X • 

z s uj  I/) 
o e r 
tn  4 o 

• X 

iij  «-*  ►- 
UJ  >«  *-* 

^ M lit  c- 
c a z z 
»-  K c 
*-  r o 
P h T 

<i  X c a 


a 

ii 

X 


« c c 
>00 


» X o 


K CVJ  ►-  C. 


131 


P^Or.ot^.  THC'SPPM(TN'''jT*''»JT''tH  -T  APrt ) 
OI'I^M^TON  IN(2J0)  ♦TOMT(l,-»1) 


c c 

ir  IT  u u 
7 2 

ir.  tr  T>  ^ 

U>  U iA  tA 
T T ir  I'i 

r;  «5  «r 

iA  tA  c-  c 

« « ut  u 


D a t. 

c c cr  c 


r c i L 
»-  *-  c c 

L ijj  r 

c c ►-«  *-• 
< « 

r r It  u 
t-  *-  C c 

; u.  u.  c u: 
' o o ft'  -» 
O c 
‘ r L c r 

. O'  ^ It  < 

> o a c/'. 
:•  e T 
■ ir  « o t-'i 

• er  c/^  s r 

• ^ ^ 

• VI  (/>  ar  a' 

I «S  < <9  < 

J -J  »-  ♦- 

I I in  1/" 

till 
till 
■■  in  in  I I 

■•  I -»  » •! 

• c c c c 

' o o ft  c‘ 


lU  t'*  < 
; CD  U'  K'  _j  1 

l Kj  1-1  O ( 

-•  I-  V c ; 

: ^ C ? M ' 

) *-<  «r  »-  e • 

O U'  O'  3 I 

• e ft'  a r > 


^ r m r c*. 

o It  3 »-  «r 

»r  r c •*  r 

•v  — ty  iri 

r z «-  » 
cvw  3 e 
# z 3 • r> 

««  • o w r 

ly  • ' 

I O ►-  Z T-  *4  1 

tLl  t-t  * tl  t-t  U’  II 

; Z'  *-t  r'  »-  D "ft  • 

: r ^ ♦ a •-  z 

■•►'•-erif  ii'Z  ►-•?“* 
.K«>’»-<iru 
r z I*  ii.  z 
:ociu.czitCo 

)o»-*^  ocr*-»c  c . 


inuT(K)5HHHSHTFT  (T»(  (JJI 
KP^KC-’ 

COMT INMt 
CONT  TN')^ 

= /■•> 


t I JdO  1 * iltvil  liu*  llltUil  Jhniulo  WtfW'JC  JO 


appendix  g 

DIGIPLT  USER'S  GUIDE 


141 


I.  INTRODUCTION 


GENERAL  DESCRIPTION 

DIGIPLT  is  an  interactive  data  curve  processor;  a collection  of  programs 
•.hat  er.ables  its  user  t.)  enter,  display,  edit,  and  operate  cn  data  ci’rv«-R. 

It  is  an  easy-to-use  software  package,  and  does  not  require  tViat  the  user 
have  extensive  knowledge  of  computer  processing. 

While  designed  priinarlly  for  use  on  the  TEKTRONIX  4014  graphics  termi- 
nal, the  user  may  take  advantage  of  many  DIGIPLT  capabilities  using  any 
interactive  computer  terminal  on-line  with  the  ASD  CDC-6600  computer  in 
INTERCOM  mode. 

The  user  should  note  that  a brief  synopsis  of  this  guide  m.ay  be  obtained 
at  the  initial  stage  of  DIGIPLT  execution  by  responding  "YES"  when  asked  by 
DIGIPLT  if  help  is  desired. 

DIGIPLT  is,  basically,  a collection  of  programs  that  are  executed  upon 
user  request.  The  request  to  execute  a particular  program  is  made  when  the 
user  specifies  one  of  the  five-character  codes  on  the  DIGIPLT  option  menu. 
Execution  of  each  of  the  programs  on  the  menu  is  completely  independent  of 
the  rest;  thus,  the  user  may  select  menu  options  in  any  order  desired. 

RESPONSE  TO  DIGIPLT  REQUESTS/QUESTIONS 

All  of  the  keyboard  entries  the  user  is  required  to  make  during  DIGIPLT 
execution  are  preceded  by  a printed  massage  (output  by  DIGIPLT)  on  the  ter- 
minal screen  specifying  the  inpun  expected.  The  response  to  DIGIPLT  requests 
Must  be  an  Integer,  or  a "YES"  or  "NO,"  depending  upon  the  nature  of  the 
request.  There  are  three  cases,  however,  when  they  must  be  real  numbers. 

1.  When  entering  axis  information  in  DGTZR. 

2.  When  entering  new  X and  Y values  when  adding  or  replacing  points 
In  EDTAL. 

3.  Wlien  entering  X and  Y values  when  creating  a data  file  in  CRTFL. 

Tf  the  request  is  for  a "YES"  or  "NO"  response,  entering  only  a "Y"  will 
n •'  -I'ffi-.e  for  1 "\'ES,"  .as  will  an  "N"  not  suffice  for  a "NO."  The  user 
in  t enter  "'.’ES"  or  "NiV " In  most  cases,  any  other  input  will  be  treated 
' f "Miy  rasoonse. 


II.  DIGIPLT  DATA  FILE  STRUCTURE 


The  purpose  of  this  section  is  to  describe  the  data  files  used  by 
DIGIPLT.  This  information  whould  be  of  use  to  all  DIGIPLT  users,  as  it  is 
essential  to  be  familiar  with  the  data  files  used  by  DIGIPLT,  since  all 
data  processed  by  DIGIPLT  must  exist  on  a DIGIPLT  data  file  in  the  correct 
format,  and  almost  all  processing  is  done  in  terms  of  these  fi'es. 

DIGIPLT  DATA  FILES 

The  user  has  available  five  separate  files  on  which  to  store  data.  Data 
is  placed  on  these  files  by  one  of  two  methods;  1)  the  user  has,  prior  to 
execution,  created  a data  file  which  is  compatible  with  the  DIGIPLT  data 
file  structure,  or  2)  the  user  creates  a data  file  through  the  use  of  the 
DIGIPLT  options  DGTZR  and/or  CRTFL.  In  the  latter  case,  the  data  file  is 
automatically  formated  in  DIGIPLT  data  file  structure.  The  user  directs 
DIGIPLT  to  use  a specific  data  file  by  entering  its  "file  number"  ("data 
file  number")  when  requested  to  do  so  by  DIGIPLT.  At  all  times  during 
execution  of  DIGIPLT,  data  files  are  referred  to  by  their  file  numbers. 

These  file  numbers  (integers)  correspond  to  actual  mass  storage  files  as 


follows : 

DIGIPLT 
FILE  NO. 

ACTUAL  MASS  STORAGE 
FILE  NAME 

1 

Tape  1 

2 

Tape  2 

3 

Tape  3 

4 

Tape  4 

5 

Tape  5 

Only  those  file^  listed  above  are  valid  DIGIPLT  data  files.  Any  attempt  to 

process  other  files  will  result  in  an  error.  Data  files  created  prior  to 

DIGIPLT  execution  should  be  on  a valid  mass  storate  file,  and  be  referred 

to  during  DIGIPLT  execution  by  its  corresponding  file  number.  Similarly,  j 

files  created  by  DIGIPLT  should  be  referred  to  by  their  actual  mass  storage 

file  name  when  used  after  completion  of  DIGIPLT  execution. 

RECORD  FORMT 

Data  exists  on  DIGIPLT  data  files  as  two  real  numbers  per  record,  in 
free-formit  (separated  by  a blank  or  a cotmna) . Each  record  corresponds  to 
one  data  point;  the  first  number  on  the  record  being  the  abscissa  value, 
and  the  second  number  the  ordinate  value  of  that  point  (X,  Y)  . 

SPECIAL  RECORDS  (FLAGS)  ! 

There  are  three  types  of  special  records  on  DIGIPLT  data  files:  1)  end- 
of-block  (E03)  flag,  also  called  end-of-curve  flags,  2)  end-of-data  (EOD)  , 

flag,  also  called  end-of-file  (EOF)  flag,  and  3)  header  block  records.  ’ 

1.  END-OF-BLOCK  (EOB)  Flag  - This  is  a record  having  values  of  { 

"9999.0"  for  both  Its  X and  Y positions.  It  is  used  to  indicate  } 
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the  end  of  a block  of  data,  where  a block  of  data  Is  treated 
by  DIGIPLT  as  one  data  (exception  - see  header  block  descrip- 
tion) curve. 

2.  END-OF-DATA  Flag  - The  last  record  in  the  data  file,  which  has 
both  an  X and  Y value  of  "8S8B.0."  It  should  be  noted  Lti.i!. 
for  the  last  block  on  the  file  there  will  be  no  £*^5  flag; 

Instead,  the  EOD  flag  (EOF  flag)  is  located  in  the  last  it-cor'f 
There  aiu..t  be  only  one  EOD  flag  on  a file;  any  data  afttr  un 
flag  Oil  DIGIPLT  data  file  will  be  ignored  and  subsequently 
destroyed.  If  the  EOD  flag  is  not  present,  an  error  will  occur 
when  attempting  to  process  the  data  file. 

3.  HEADER  BLOCK  RECORDS  - The  first  block  on  all  DIGIPLT  data 
files  contains  information  on  the  contents  of  the  file  and  is 
used  by  DIGIPLT  for  scaling  purposes.  This  block  is  three 
records  (points)  in  length,  and  contains  the  following; 

REG  X POSITION  Y POSITION 

1 XMAX  YMAX 

2 XMIN  YMIN 

3 EOB 

where:  XMAX  and  YMAX  are  the  maximum  X and  Y values,  respec- 
tively, of  the  data  points  on  the  file;  XMIN  and  YMIN  are  the 
minimum  X and  Y values,  respectively,  of  the  data  points  on  the 
file. 

It  should  be  noted  that,  due  to  the  header  block,  the  first  actual  data 
curve  on  a file  is  actually  the  second  data  block.  The  user  should  keep 
this  in  mind  when  DIGIPLT  asks  for  "BLOCK  NO."  and  "CURVE  NO."  from  the 
ncer  - a data  curve's  "BLOCK  NO."  will  always  be  one  greater  that  its 
'CURVE  NO." 

PERMANENT  FILES 

If  the  user  desires  to  use  data  located  on  a permanent  file  cataloged 
,rior  to  DIGIPLT  execution,  he  must  copy  the  file  to  one  of  the  DIGIPLT 
iota  files  before  execution  of  DIGIPLT.  If  a permanent  file  is  attached  as 
i-'e  of  the  OTGrPLT  dat."!  files,  an  unrecoverable  error  will  occur  in  DIGIPLT 
execution  due  to  the  fact  that  DIGIPLT  writes  on  the  data  files  it  processes, 
;hlch  IS  an  illegal  operation  on  a permanent  file.  The  user  m.ay,  hoa^ever, 
'.-qu'.  f. t i a UjUJPLT  data  file  be  assigned  to  a permanent  file  device  and 
...c  — .'.'.if  il..;  altec  iKCTi’LT  execution. 


I 

m,  MENU  OPTIONS  ; 

I 

The  following  section  elaborates  on  each  of  the  user-selectable  programs  ; 

of  DIGIPLT  (menu  options).  A menu  option  is  executed  when  the  user  enters  ! 

its  5-character  option  code  in  response  to  the  DIGIPLT  query  "ENTER  5-CHAR. 

OPTION  CODE." 


ALPHA-NUMERIC  DATA  CURVE  DISPLAY  (DSPAL) 

The  function  of  DSPAL  is  to  produce  a listing  of  the  abscissa  and  ordi- 
nate values  of  points  on  specific  data  blocks.  This  display  will  indicate 
tho  block  t.u/nbei , f-oif.i.  number,  and  X and  Y values  for  all  bloik  ■ sp  ified 
to  be  displayed  on  a given  data  file. 

Since  DSPAI.  displays  are  in  terms  of  blocks,  the  user  may  displa/  tlie 
header  block  of  each  file.  Therefore,  the  first  actual  data  cur're  Oi  the 
fil«‘  will  be  block  no.  2.  When  asked  to  input  the  block  numbers  to  be 
displayed,  an  entry  of  "0"  for  the  first  block  number  will  cause  all  blocks 
on  the  specified  data  file  to  be  displayed.  If  the  user  chooses  to  display 
only  certain  blocks,  they  should  be  entered  in  ascending  order  by  their 
block  numbers.  A "0"  must  be  entered  as  the  last  block  to  be  displayed. 


TYPE  IN  5-CHAP.  FUHCTIOI 
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GRAPHIC  DATA  FILE  DISPLAY  (DSPGR) 


The  function  of  the  DSPGR  option  is  to  produce  a plot  on  the  terminal 
screen  of  specified  data  curves  on  a specific  data  file.  The  plot  produced 
Is  a simply  X-Y  graph  of  the  data  curves,  and  enables  the  user  to  examine 
the  data  in  graphic  form. 

Since  DSFAL  displays  are  in  terms  of  data  curves,  a "1,"  for  example, 
entered  as  the  data  curve  to  be  displayed  will  result  in  a plot  of  the  first 
actual  data  curve  on  the  specified  file.  As  in  DSPAL,  a "0"  e.-?tered  for  the 
first  data  curve  to  be  displayed  will  result  in  the  display  of  all  data 
curves  on  the  file.  A "0"  must  be  entered  as  the  last  curve  to  be  displayed. 

The  criteria  for  scaling  of  the  displayed  graph  is  determined  by  the 
method  in  which  the  user  chooses  to  display  the  data  curves.  ilf  he  elects  to 
display  all  curves  on  a file  by  entering  a "0"  as  the  first  curve  to  be  dis- 
played, scaling  of  the  graph  will  be  in  terms  of  the  information  in  the 
header  block  of  the  data  file.  If  a specific  curve(s)  to  be  displayed  is 

entered,  scaling  of  the  graph  v;ill  take  place  in  relation  to  the  minimum 

and  maximum  abscissa  and  ordinate  values  of  the  data  to  be  displayed. 

The  user  has  the  option  of  displaying  the  data  curves  on  each  plot  by 
two  methods: 

1.  Consecutive  points  on  each  curve  connected  by  straight  lines. 

2.  Consecutive  points  on  each  curve  connected  by  straight  lines, 

with  a symbol  drawn  at  each  point.  (A  different  symbol  is  used 
for  each  individual  curve  on  a given, plot.) 

The  user  selects  the  method  by  responding  "NO"  or  "YES,"  corresponding  to 
(1)  and  (2)  respectively,  when  asked  if  symbols  are  to  be  displayed. 
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ALPHA-NUMERIC  DATA  FI1.E  EDITOR  (EDTAL) 

The  purpose  of  EDTAL  is  Co  provide  Che  user  wich  Che  capabllicy  of 
edicing  daCa  blocks.  As  in  DSPAL,  EDTAL  refers  Co  Che  daCa  file  in  Cerms 
of  blocks,  Chus  allowing  Che  edicing  of  Che  header  block  on  a file. 

The  user  has  Chree  edicing  funcCions  available;  1)  addicion  of  new 
poinCs  Co,  2)  deleCion  of  poinCs  from,  and  3)  replacing  of  poinCs  wiCh  new 
ones  on  a daCa  block.  These  funcCions  are  effecCed  by  enCering  a one-leCCer 
edic  command  when  asked  by  DSPAL  Co  do  so.  Valid  edic  commands  are  as 
follows: 

1-LeCCer  Command  FuncCion 

A Add  poinC(s)  Co  specified  block 

D DeleCe  poinC(s)  from  specified  block 

R Replace  poinC(s)  on  a specified  block 

The  user  is  asked  Co  enCer  the  poinC  number  to  be  processed  by  EDTAL. 
This  poinC  number  enCered  here  should  be  Che  poinC  number  on  Che  display 
corresponding  Co  che  posicion  the  user  wishes  Co  begin  edicing,.  In  Che 
case  of  deleting  or  replacing  points,  the  point  specified  will  actually  be 
affected  (example:  point  no.  specified  is  5,  edit  command  was  D,  point 
no.  5 will  be  deleted).  In  the  addition  of  point(s),  the  new  point(s)  will 
be  inserted  in  Che  block  immediately  after  the  point  number  specified. 

Both  before  and  after  executing  an  edit  command,  the  user  has  the  option 
of  displaying  the  block  he  wishes  to  edit.  All  "point  numbers"  referred  to 
in  this  discussion  are  those  which  appear  on  this  display.  (This  display 
is  of  the  same  nature  as  that  produced  by  DSPAL.) 

When  adding  or  replacing  points,  the  user  must  enter  the  new  X and  T 
values  for  each  point.  They  must  be  entered  as  real  numbers,  separated  by 
a blank  or  a comma. 

The  user,  when  specifying  a deletion,  has  the  option  of  deleting  the 
entire  block.  If  a "YES"  response  is  given  when  the  user  asks  if  he  wishes 
this  done,  the  entire  block  will  disappear  from  the  data  file.  Note  that 
if  the  last  block  of  the  file  is  deleted,  the  last  record  (end-of-block 
flag)  of  the  preceding  block  must  be  changed  from  "9999.0  . 9999.0"  to 
"8888.0  . 8888.0"  Co  assure  normal  processing  of  that  file  and  prevent  the 
occurance  of  errors.  In  the  same  respect,  the  user  may  add  an  entire  block 
to  the  file  by  adding  points  after  the  end  of  block  or  end  of  data  (file) 
flag,  specifying  as  the  last  X-Y  pair  the  end-of-block  or  end-of-data  flag 
accordingly.  In  determining  proper  use  of  the  above  flags,  it  is  important 
to  remember  that  any  data  appearing  after  an  end-of-data  flag  ("8888.0  . 
8888.0")  will  be  ignored  and  consequently  deleted  from  Che  data  file. 
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OR  A CONHA.  FOR  the  POINT  NO.S  SPECIFIED 


DATA  CORVE  INTERPOLATION  (INTRP) 


The  function  of  the  INTRP  option  is  to  interpolate  points  between  the 
existing  points  on  a specified  data  curve(s).  As  in  DSFAL  and  DSPGR,  a 
"0"  entered  for  the  first  curve  to  be  processed  will  result  in  interpolation 
on  all  curves  on  the  specified  data  file;  if  specific  curve(s)  are  to  be 
processed,  a "0"  must  be  entered  as  the  last  curve  to  be  interpilated. 

The  user  is  asked  to  enter  the  number  of  intervals  to  be  interpolated 
between  the  original  points  on  the  data  curve(s).  In  order  for  Interpola- 
tion to  occur,  this  number  must  be  in  the  range  of  2 - 10  inclusive.  INTRP 
will  interpolate  and  insert  between  the  original  points  of  the  data  curve(s) 
the  number  of  points  equal  to  one  less  than  the  number  of  intervals  speci- 
fied. (Thus,  if  an  interval  of  "1"  was  specified,  no  points  would  be 
interpolated.) 

Due  to  limited  storage  space  in  the  computer,  individual  data  curves 
are  limited  to  having  a maximum  of  900  points  each.  Should  the  number 
of  intervals  specified  result  in  a curve  of  length  greater  than  900  points, 
the  largest  interval  which  will  contain  the  curve  to  the  maximum  903  points 
is  automatically  used  by  INTRP. 

In  the  actual  interpolation,  INTRP  uses  a 3d  degree,  piece-wise  fitting 
of  the  original  data  points  to  determine  the  value  of  interpolated  points. 
The  abscissa  values  of  the  interpolated  points  are  determined  by  the  number 
of  intervals  specified  and  the  abscissa  values  of  the  original  points;  their 
ordinate  values  are  determined  by  the  abscissa  values  and  the  fit.  It  is 
mandatory  that  any  curve  being  interpolated  by  continuous,  and  the  abscissa 
values  be  Increasing  as  their  point  numbers  increase. 
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DATA  FILE  STATISTICS  AND  COPyiNG  (FILES) 


The  function  of  the  FILES  option  Is  to  provide  the  user  with  the  status 
of  each  of  the  five  data  files  at  his  disposal,  and  enable  the  copying  of 
one  data  file  to  another. 

The  number  of  data  curves,  the  number  of  points  on  each  curt/e,  and  the 
total  number  of  points  contained  in  the  file  are  given  for  each  file.  If  a 
file  contains  no  data,  the  message  "FILE  EMPTY"  is  given  for  that  file. 

In  copying  files,  any  valid  data  file  may  be  copied  to  any  other  valid 
data  file.  If  the  user  attempts  any  copy  involving  non-existent  files,  an 
error  message  is  printed  and  no  copying  occurs.  If  a file  that  is  specified 
to  be  copied  to  contains  data  prior  to  the  copy,  a warning  is  given;  and  if 
the  user  elects  to  proceed  with  the  copy,  any.  data  previously  on  the  file 
will  be  destroyed. 


T'.i’t  In  S-CH«W.  FUNCTION  lAME 


r 

DATA  FILE  CREATION  FROM  KEYBOARD  ENTRY  (CRTFL^ 

The  purpose  of  the  CRTFL  option  is  to  provide  the  user  with  the  capability 
of  creating  data  curve(s)  by  entering  X-Y  values  from  the  terminal  keyboard. 
Data  curves  may  be  created  on  any  valid  DICIPLT  data  file.  (A  warning  will 
be  Issued  If  an  attempt  Is  made  to  create  data  curves  on  a file  currently 
containing  data.) 

The  process  of  the  creation  is  a point -by-polnt  entry  of  data  curves. 

An  end  of  block  flag  entered  (9999.0  . 9999.0)  as  the  X and  Y values  results 
in  the  beginning  of  a new  curve  on  the  file.  The  last  point  entered  must  be 
an  end-of-data  flag  ("8888.0  . 8888.0");  this  results  In  the  completion  of 
the  data  file  creation.  CRTFL  will  continue  to  ask  for  X-Y  values  until 
the  end-of-data  flag  Is  entered.  An  end-of-block  flag  must  not  be  entered 
on  the  last  block  to  be  created.  Each  point  must  be  entered  as  an  X and 
a Y value,  separated  by  a blank  or  cooma,  on  the  same  line. 

The  header  block  required  on  all  DIGIPLT  data  files  Is  automatically 
created  and  placed  in  the  correct  location  on  the  file  created. 


I 
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This  option  should  be  executed  only  if  using  a TEKTRONIX  4014  graphics 
terminal,  with  a TEKTRONIX  4954  digitizing  tablet.  Any  mention  of  the 
"digitizer"  in  the  following  discussion  refer  to  the  4954  tablet  and  con- 
troller. 

TABLET  OPERATION 

This  option  provides  the  user  with  the  capability  of  entering  data  curves 
by  digitizing  graphs.  If  it  is  intended  to  use  this  option,  it  is  advisable 
to  turn  on  the  tablet  controller  before  beginning  execution  of  DIGIPLT. 

Secure  the  graph  to  be  digitized  to  the  tablet  with  four  small  pieces 
of  adhexive  tape.  It  is  advisable  to  place  the  graph  near  to  the  center  of 
the  tablet,  as  the  usable  area  of  the  tablet  does  not  necessarily  extend  to 
its  edges.  DIGIPLT  automatically  corrects  for  skew  of  the  graph  with  respect 
to  the  tablet.  While  digitizing  do  not  lift  the  mouse  from  the  tablet.  If 
this  happens  while  the  digitizer  is  armed  (ready  light  on  controller  is  on), 
the  ready  status  of  the  digitizer  will  be  lost.  In  this  case,  return  the 
mouse  to  the  tablet  and  move  it  around  until  the  ready  status  returns.  Do 
not  key  the  mouse  unless  the  tablet  is  armed.  To  digitize  a point,  place 
the  crosshairs  of  the  mouse  over  the  point,  and  key  the  mouse  once.  Do  not 
move  the  mouse  while  it  is  keyed.  DIGIPLT  automatically  arms  the  tablet 
when  it  expects  input  from  the  digitizer. 

GRAPH  FORMAT 

Any  graph  that  is  to  be  digitized  must  have  an  X (abscissa)  axis,  and 
a Y (ordinate)  axis  perpendicular  with  respect  to  the  X axis,  associated 
with  the  graph  (i.e.  must  have  a cartesian  coordinate  system). 

It  is  also  required  that  there  be  associated  with  these  axes: 

1.  an  abscissa  and  ordinate  value  at  the  origin. 

2.  an  abscissa  value  at  the  maximum  extent  (pos.  sense)  of  the  X 
axis. 

3.  an  ordinate  value  at  the  maximum  extent  (pos.  sense)  of  the  Y 
axis. 

These  values  may  be  arbitrary.  The  points  digitized  will  have  coordinates 
valued  with  respect  to  the  above.  To  establish  the  relation  between  the 
values  of  the  points  above  to  the  coordinate  values  output  by  the  tablet, 
DGTZR  will  ask  the  user  to  type  in  the  values,  (1,  2,  and  3 aoove)  and  then 
digitize  the  corresponding  points  on  the  graph. 

At  this  point,  DGTZR  will  ask  the  user  to  input  (digitize)  the  curve 
with  the  mouse,  entering  a "D"  from  the  keyboard  when  finished  digitizing 
that  curve.  When  entering  a curve,  move  the  mouse  along  the  desired  curve, 
keying  it  at  the  locations  where  actual  points  are  to  be  t 'ken.  Begin  at 
the  left-mcst  point  on  the  curve  and  proceed  to  its  end.  Ine  output  of  the 
digitizer  is  echoed  on  the  terminal  screen  (appears  as  "strange"  character 
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strings).  If  miny  points  are  being  digitized,  be  careful  that  the  terminal 
screen  does  not  become  "full,"  causing  the  terminal  to  stop  receiving  data 
from  the  digitizer.  When  the  curve  has  been  digitized,  type  a "D"  and  a 
carriage  return  on  the  keyboard.  If  at  this  time  there  Is  another  curve  to 
be  digitized,  respond  "YES"  when  asked  by  DGTZR  if  another  curve  is  to  be 
A "I'.O"  response  here  will  result  in  final  processing  of  the  da>:a 
entered  and  exit  from  DGTZR. 


EXIT  FROM  MENU  OPTION  SELECTION  MODE  (FINSH) 


The  purpose  of  the  FINSH  option  is  to  exit  from  the  user-selectable 
option  mode.  This  option  is  executed  when  the  user  has  completed  processing 
all  data.  Upon  exiting,  DIGIPLT  will  ask  the  user  if  he  is  finished.  A 
"YES”  reply  given  at  this  point  will  result  in  the  termination  of  DIGIPLT 
execution.  All  data  will  remain  on  the  DIGIPLT  data  files,  for  the  user  to 
do  with  as  he  wishes  (example:  catalog  as  permanent  file,  dispose  to  a 
batch  terminal,  or  use  as  input  to  another  program).  A "NO"  response  will 
result  in  the  return  of  the  user-selectable  option  mode,  and  the  user  may 
continue  processing. 
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IV.  USE  OF  DIGIPLT  FROM  TCRMItlAL  OTHER  THAN  TEKTRONIX  4014 

DIGIPLT,  while  primarily  designed  to  be  executed  using  a TEKTRONIX 

executed  using  any  terminal  on-line  with  the  CDC-660D 
(INTpRCOM).  Two  restrictions  are  observed; 

1.  DO  not  SELECT  DGTZR  OPTION. 

2.  DO  NOT  SELECT  DSPGR  OPTION. 

Ihc  procedures  for  use  of  all  DIGIPLT  options  other  than  chose  above  from 
a termtnal  other  than  the  TEKTRONIX  4014  are  exactly  as  outlined  In  this 
guide. 
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V. 


GENERAL  INFORMATION 


EXECUTION  OF  DIGIPLT 

The  DIGIPLT  program  is  maintained  in  absolute  binary  form  on  the  AFML 
program  library.  To  execute  DIGIPLT,  the  user  must  login  to  INTERCOM  on 
the  rDC-6600,  attach  DIGIPLT  from  the  library,  and  execute.  Since  DIGIPLT 
is  In  absolute  form  it  requires  no  additional  supporting  software. 

DIGIPLT  is  attached  and  executed  by  the  following  commands: 

ATTACH, DIGIPLT, DIGIPLT , ID=M7 54321 . 

DIGIPLT. 

These  commands  must  be  given  exactly  as  shown.  Any  variations  will  result 
in  execution  errors. 

The  user  should  become,  if  not  already,  familiar  with  the  use  of  CDC-6600 
INTERCOM.  Information  on  its  use  is  available  from  AFML/DOC. 

IF  DIGIPLT  BOMBS 


Should  the  user  cause  DIGIPLT  to  bomb,  all  data  files  are  preserved  in 
the  state  they  were  in  immediately  before  the  termination  of  execution. 
INTERCOM  commands  may  be  used  to  examine  data  (ex.:  COPYSBF,  TAPEl,  OUTPUT) 
attempt  to  ascertain  the  cause  of  the  error,  if  not  known  already.  The  user 
may  elect  to  use  DIGIPLT  to  examine  the  data.  In  this  case,  it  is  strongly 
advisable  to  select  the  FILES  option  before  any  other,  as  this  will  give  the 
user  the  most  information  on  the  data  files. 

REPORTING  BUGS  AND  ERRORS 


Should  the  user  discover  a bug  or  error  in  DIGIPLT,  he  is  asked  to  docu- 
ment and  report  the  problem  to  AFML/DOC.  The  documentation  need  only  be  the 
terminal  output  sufficient  to  show  the  existence  of  the  bug  or  error,  or  the 
equivalent  explanation. 

ADDITIONAL  INFORMATION 

Manuals  on  the  use  of  INTERCOM,  the  TEKTRONIX  4014  terminal,  and  TEKTRONIX 
4954  digitizing  tablet  are  available  at  AFML/DOC,  should  the  user  desire  more 
information  on  these. 


VI.  SPECIFICATIONS 


APPENDIX  H 

GENERAL  USER  PROGRAMS  FOR  THE  CDC  6600 
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